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O ■ ABSTRACT 

(N; 

Aims. We describe our newly developed approach to detailed abundance analysis from integrated-light high-dispersion spectra of 
■ star clusters. As a pilot project, we measure abundances of several Fe-peak, a- and neutron capture elements from spectra of three 

" globular clusters (GCs) in the Fornax dwarf spheroidal galaxy, obtained with UVES on the ESO Very Large Telescope. 

Methods. We divide the cluster colour-magnitude diagrams into about 100 bins and use the Kurucz ATLAS9 and SYNTHE codes to 
compute synthetic spectra for each bin. Stellar parameters are derived empirically from Hubble Space Telescope data for the brighter 
stars, while theoretical isochrones are used for extrapolation below the detection limit. The individual model spectra are co-added and 

Othe abundances are iteratively adjusted until the best match to the observed spectra is achieved. 
Results. We find [Fe/H] = -2.3, -1.4 and -2.1 for Fomax 3, 4 and 5, with estimated ±0.1 dex uncertainties. Fornax 3 and 5 are thus 
similar in metallicity to the most metal-poor Milky Way GCs and fall near the extreme metal-poor end of the field star metallicity 
distribution in Fornax. The [a/Fe] ratios, as traced by Ca and Ti, are enhanced with respect to the Solar composition at the level of 
i-C 1 ~ +0.25 dex for Fornax 3 and 5, and possibly slightly less (~ +0.12 dex) for Fornax 4. For all three clusters the [Mg/Fe] ratio is 

significantly less elevated than [Ca/Fe] and [Ti/Fe], possibly an effect of the abundance anomalies that are well-known in Galactic 
JL , GCs. We thus confirm that Mg may be a poor proxy for the overall a-element abundances for GCs. The abundance patterns of heavy 

elements (Y, Ba and Eu) indicate a dominant contribution to nucleosynthesis from the r-process in all three clusters, with a mean 

+jf ■ [Ba/Eu] 0.7, suggesting rapid formation of the GCs. 

Conclusions. Combining our results with literature data for Fornax 1 and 2, it is now clear that four of the five Fornax GCs fall in the 
range -2.5 < [Fe/H] < -2, while Fornax 4 is unambiguously and substantially more metal-rich than the others. The indications that 
abundance anomalies are detectable in integrated light are encouraging, particularly for the prospects of detecting such anomalies in 
£N| ' young, massive star clusters of which few are close enough for individual stars to be observed in detail. 
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q^! 1. Introduction type la SNe dTinslevl[T979t iMatteucci & Greggid fl986h . The 

ancient stellar populations in the Milky Way halo and bulge, as 

^ ■ The chemical composition of stellar populations as a function well as giant e arly-type galaxies, are typically enhanced in the 

of time and location holds vital clues to the evolution of their ^-eleme nts, indicating t hat star form ation was a relatively rapi d 

host galaxies. The distribution of metalhcities for G-dwarfs process dWorthevlll998l:lTrager et al.ll2000l:lThomas et alT2 005). 

in the Galactic disc implies signincajrtgas infall over the Further in f ormation about time sca i es and nucleosynthetic 

. lifetime of the disc (van den Bergh 1962; Pagel & Patchettl si tes comes from the abundances of neutron-capture elements 

r £ -119751; | Kotonevaetal. | |2002j; | Naab & Ostrikeij |2006j), and (McWillii^Clll) and even within the Fe-peak group, different 

& ; age-metalhcity and Galactocentnc distance-metalhcity trends elements display significant abundance variations that may 

provide additional important constraints on modelsfor hold important clues to the detailed nucleosynthetic histories of 

Galactic chemical e volutio n (jMatteucci & Francois 1989; stellar populations 



Pagel & Tautvaisiene 1995; IChiappini et all I1997L I200U 



Marcon-Uchida et al.1 12010I) . An even more detailed picture Although a vast body of work exists on the chemical com- 



comes from the abundance patterns of individual elements, position and evolution of both Galactic and extragalactic stellar 

which are sensitive to the time scales of star formation. Perhaps populations, there is a striking difference between the level of 

the best-known example is the ratio of a-capture to Fe-peak detail that can typically be achieved in these two domains. In 

elements, which is an indicator of the relative contributions to the Milky Way, high-dispersion spectroscopy of individual stars 

nucleosynthesis by (short-lived) core-collapse supernovae vs. provides a highly d etailed picture of the abundance patterns o f 

individual elements (E dvardsson et al. 1993; Re ddv et al.l 2006). 

Send offprint requests to: S. S. Larsen, e-mail: either via measurements of the equivalent widths of individual 

s.larsen@astro.ru.nl lines or, when line blending gets too severe, spectral synthesis. 

* Based on observations made with ESO Telescopes at the La Silla Such measurements can be made for stars of essentially all ages, 

Paranal Observatory under programme ID 078.B-0631(A) although the most detailed measurements remain constrained to 
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the relative vicinity of the Sun (though recently pu shed to dwarf- 
andsubgiant stars in the bulge via microlensing; Bensbv et al. 
1201 lb . For external galaxies, constraints on chemical compo- 
sition must typically rely either on measurements of emission 
lines in H n regions (for star-forming galaxies) or on relatively 
low-dispersion spectroscopy of integrated stellar light. In most 
cases, interstellar gas abundances are only available for a limited 
number of elements (mainly O, N, S) and rely on em pirical and 
uncer tain calibrations of line strength vs. abundance ( Sta sinskal 
l2010l) . Furthermore, emission lines only provide a snapshot of 
the present-day gas chemistry a nd may be biased by e ffects such 
as segregation onto dust grains dPeimbert & Peimbert 2010). 

Integrated-light measurements of stellar absorption features 
are typically carried out by means of pre-defined sets of 
abso rption-line indices, t he best-known being the Lick/IDS sys- 
tem (IWorthev et al.1 H994I) . These indices are relatively broad 
(~ 10 - 50A) and well-suited for measuring strong absorption 
features in the integrated light of early-type galaxies, where the 
spectral resolution is in any case limited by the internal velocity 
dispersions of the galaxies. Although the various indices have 
different sensitivities to different elements (such as Fe, Mg and 
the hydrogen Balmer lines), none of them cleanly measures in- 
dividual absorption lines. Nevertheless, the Lick/IDS system has 
been very successful in providing constraints on (luminosity- 
weighted mean) ages and overall metallicities of extragalactic 
stellar populations, and e ven on some abundance patterns such 
as overall f g/Fe] ratios (iKuntschner j UOQOt iTrager et aD [2000: 
iPuzia et al.|[2005t Thom as et al. lool . Physical properties such 
as ages and metallicities are normally inferred from such data 
by comparison with predictions by "simple stellar population" 
(SSP) models, either by simply plotting suitable combinations 
of index measurement s on top of a model grid or in a more 
automated fashion (e.g.lCenarro et al.ll2007 : Graves & Schiavonl 
2008; iProctor et al.l l2004t IStrader & Brodid 120041) . In some 
cases, more sophisticated inversion techniques have also been 
appli ed to infer star formation histories from spe c troscopic data 
dHeayens et al.ll2004t IcTd Fernandes eTaTl 120051: IWalcher et all 
120101) . 

Unlike the spectra of giant elliptical galaxies, the spectra of 
most star clusters are velocity broadened by only a few km s _1 , 
allowing observations at much higher spectral resolution. With 
efficient high-dispersion spectrographs on 8-10 m class tele- 
scopes, such as VLT/UVES and Keck/HIRES, high S/N spectra 
can be obtained for star clusters well beyond the Local Group. 
Such observations have the potential to provide a vastly more 
detailed picture of chemical abundance patterns compared to 
techniques developed for integrated-light spectra of galaxies. 
Furthermore, star clusters provide discrete sampling points in 
time and space, even for stellar populations that contribute only 
a small fraction of the integrated light (such as metal-poor stellar 
halos). However, the tools required to interpret such data have, 
until recently, not been available. 

Over the past decade, our group has been active in collect- 
ing high-dispersion spectra of extragalactic star clusters. These 
have been used primarily for dynamical mass measurements 
(Larsen et al.ll2002tlLarsen & Richtledl2004t ILarsen et al.ll2004l 
Stra der et al.ll2009l.l201 lb . but in many cases the S/N is also ad- 
equate for abundance analysis. With the next generation of ex- 
tremely large telescopes, such data will be within reach at the 
distance of the Virgo cluster and beyond. Time is therefore ripe 
to develop the techniques required to extract the information 
available in such data. 

It should be emphasized that the use of star clusters as probes 
of stellar populations does come with its own set of caveats. 



The relative numbers of globular clusters and field star s vary 
strongly, both from galaxy to galaxy and within galax ies (iHarrisI 
I199U iHarris & Harris! Ilool iBrodie & Stradeiil2006h . The rea- 
sons for these variations are poorly understood, but are likely 
due to a combination of dynamical evolutionary effects and vari- 
ations in the formation efficiency of massive star clusters rela- 
tive to field stars. While clusters may be used to identify bursts 
of star formation, e.g., due to major mergers, inferring the in- 
tensities of such bursts is therefore far less straight forward. In 
the context of chemical abundance analysis, it is also important 
to note that some elements are known to exhibit peculiar abun- 
dance patterns in globular clusters. These anomalies mostly af- 
fect the light elements He, C, N, O, Na, Mg and Al, and sug- 
gest t hat self-enrichmen t has taken place within globular clus- 
ters dGratton et alj|2012h . However, apart from these specific el- 
ements, the abundance patterns of globular clusters are generally 
similar to those seen in field stars, and at least in the Galactic 
disc both open clusters and field stars paint a c onsistent pi cture 
of age-metallicity relations and radial gradients dFrielf 1995). 

Compared to high-dispersion spectroscopy of individual 
stars, detailed abundance analysis from the integrated spectra of 
star clusters introduces a number of additional challenges: First, 
the integrated optical light generally has significant contributions 
from stars across the Hertzsprung-Russell (H-R) diagram. In 
some cases this complication can be circumvented by choosing 
the wavelength range carefully. In particular, the near-infrared 
spectra of young (:$ 100 Myr) star clusters are dominated by red 
supergiant stars, and for such objects it may suffice to model the 
IR spectr um of an entire cluster as that of an individu al, red su- 
pergiant dLarsen et al. 1 120061 [20081: iDavies et al.ll20loh . Second, 
although the features in spectra of star clusters are much less 
broadened by internal velocity dispersions compared to, say, gi- 
ant elliptical galaxies, the broadening can still lead to significant 
blending of spectral features. Measurements of the equivalent 
widths of individual lines therefore become more problematic 
than for single stars, requiring an increased reliance upon spec- 
tral synthesis. In recent years, there has been a general move 
towards computing SSP models at increasing spectral resolu- 
tion and wider w avelength coverage, and grids of such mod- 
els ar e available dGonzalez Delgado et al. 2005; Coelho "et al.l 
|2007|) . However, for detailed abundance analysis it becomes im- 
practical to rely on such pre-calculated grids, which would have 
to be tabulated for a large range of ages, overall metallicities, 
horizontal-branch morphologies, detailed chemical abundances, 
etc., and at a sufficiently fine spacing to allow constraints on the 
abundance patterns at the ~ 0. 1 dex level or better. 

Our aim has been to develop a general method that is not 
so specialized that it is just applicable to high-dispersion, high 
S/N data of simple stellar populations with relatively modest 
internal velocity dispersions. In future applications, we may 
also apply it to less ideal cases, such as data from spectro- 
graphs that deliver somewhat lower resolution (but are highly 
efficient, e.g., X-Shooter on the ESO VLT). Conceptually, our 
ap proach has several aspect s in c ommon with that introduced 
bv iMcWilliam & Bernsteinl d2008l) and further developed by 
IColucci et alJ (120091 1201 ll 120121) . In both cases, the integrated 
light measurements are interpreted based on the detailed H-R 
diagrams of the clusters, combined with model atmosphere cal- 
culations and synthetic spectra of the kind that have traditionally 
been used for abundance analysis of individual stars. However, 
we model larger sections of the integrated spectra, solving si- 
multaneously for (in principle) arbitrary combinations of indi- 
vidual element abundances. This allows us to use information 
contained in many weak lines that may not be measurable indi- 



2 



S. S. Larsen et al.: High-dispersion spectroscopy of the Fornax dSph GCs 



vidually, but still constrain the abundance measurements. A fur- 
ther motivation is that we may also, in later applications, apply 
this type of analysis to medium-resolution spectra that will be 
more strongly affected by blending, especially for higher metal- 
licities. This will clearly involve trade-offs in the level of detailed 
information that can be extracted from the data, but we believe 
that the detailed spectral synthesis approach still has the poten- 
tial to extract much more detailed information compared to, e.g., 
Lick/IDS index measurements. 

In this paper we introduce our approach and apply it to three 
of the globular clusters in the Fornax dwarf spheroidal galaxy. 
We start by briefly discussing previous work on the metallic- 
ities of the Fornax GCs (Sect. [2]). Our spectroscopic observa- 
tions, their basic reduction, and the derivation of physical pa- 
rameters for stars in the clusters are then discussed in Sect. [3] 
followed by a detailed description of the procedure by which we 
derive abundances from integrated-light spectroscopy (Sect. @J. 
In Sect. [5] we present our results and carry out a number of tests 
to assess the various uncertainties. In Sect.[6]we discuss our mea- 
surements in the context of other abundance determinations for 
metal-poor stars in dwarf galaxies and the Milky Way, and our 
main conclusions are summarized in Sect. [7] 



2. Globular clusters in the Fornax dwarf spheroidal 
as a test case 

The globular clusters Fornax 3, Fornax 4 and Fornax 5 in the 
Fornax dwarf spheroidal galaxy (Hod ge|[l96ll) were chosen as 
an ideal first test sample for several re asons: 1) With core radii o f 
2"-3" and half-light radii of 5"-10" dMackev & Gilmorel2003l) . 
they are sufficiently compact that integrated-light spectra can be 
obtained relatively easily. 2) At the same time, they are suffi- 
ciently nearby that resolved photometry from the Hubble Space 
Telescope ( HST), reaching below the main sequence turn-off, 
is available dBuonanno et al. 1998, 119991) . Finally, although this 
was not known at the time our observations were planned, de- 
tailed abundance analysis has been carried out for individual 
red giant branch (R GB) stars in one of the clusters (Fornax 3; 
iLetarte et alj 120061) . providing a welcome comparison for our 
integrated-light results. 

Previous determinations of the metallicities of the Fornax 
GCs have b een extensively reviewed and discussed by 
other authors dBuonanno et al.lll998lll999t IStrader et all [2003 ; 
ILetarte et all 120061: iMottini & Wallersteinl 120081) . While most 
studies agree that Fornax 1, 2, 3 and 5 are very metal-poor 

with [Fe/H] 2, there are still substantial uncertainties. Past 

studies have relied on a variety of techniques, including inte- 
grated colours, resolved photometry and low- or moderate dis- 
persion integrated spectroscopy, most of which require interme- 
diate calibration steps that are not always well established at 
such low metallicities. From measurements of Lick/IDS indices 
on Keck/LRIS spectra, IStrader et alj {2003) derived somewhat 
higher metallicities of [Fe/ H] ~ —1.8 f o r Forn ax 2, 3 and 5. 
The UVES spectroscopy of ILetarte etal] d2006l hereafter L06) 
instead yielded very low metallicities of [Fe/H] = -2.5, -2.1 
and -2.4 (with estimated +0.1 dex errors) for a few individual 
red giants in Fornax 1, 2 and 3, respectively. 

The case of Fornax 4 is somewhat more controversial. 
Generally, this cluster has been thought to be more metal-rich 

than the other four, with typical values of [Fe/H] 1.4 being 

derived from integrated photo metry and/or low-dispersion spec- 
troscopy (see references in iBuonanno et alj 119991 : IStrader et alj 
120031) . IStrader et alJ d2003l) found [Fe/H] = -1.5 for Fornax 4, 




Fig. 1. Image of Fornax 4 from the UVES slit- viewing camera at 
the start of an East- West scan. The 10" long UVES slit is visible 
at its initial position, one half-light radius (5'.'5) from the centre 
of the cluster. 



and argued (given their relatively high metallicity estimates for 
the other clusters) that the difference between Fornax 4 and the 
other GCs might be smaller than suggested b y prev ious stud- 
ies. From HST photometry, Buonann o et al.1 ([T999) also con- 
cluded that the metallicity of Fornax 4 is not significantly higher 
than that of the other clusters. In their case, this conclusion was 
based on their finding that the RGB slope was similar to that 
of the other Fornax GCs, suggesting a similar (low) metallicity 
of [Fe/H] ~ -2. They did note that this was in disagreement 
with spectroscopic observations available at the time, and that 
the red horizontal branch of Fornax 4 would also be atypical of 
such a metal-poor cluster although this might be in part due to a 
younger age. 

Finally, we note that Fornax 5 is known to host a planetary 
nebula that w as discovered in the same data used in this paper 
dLarsenll2008l) . The PN is peculiar in having no detectable H 
Balmer emission lines, in spite of strong [O m] emission, and is 
one of only two known PNe in the Fornax dSph. 

3. Observations and data analysis 

3.1. UVES spectroscopy 

We obtained integrated-l ight spectra of Forn ax 3, 4 and 5 with 
the UVES spectrograph dPekker et alJl2000h on the ESO Very 
Large Telescope on 19 Nov and 20 Nov, 2006. The observations 
were made with the UVES red arm, using cross-disperser #3 and 
a slit width of 1". This set-up yielded a resolving power R ~ 
40 000 and a spectral range of 4200 A - 6200 A, with a gap 
from 5170 A - 5230 A between the two CCD detectors. The 
seeing varied between (X'5 and 1"5 during the two nights, but the 
actual seeing was of little importance for our purpose. Weather 
conditions were mostly clear, with some clouds near the horizon 
that did not significantly affect our observations. 

To sample the integrated light, we scanned the UVES slit 
across each cluster in both the North-South and East- West di- 
rections. Each scan had an exposure time of 2400 s, and during 
this time the slit was moved across the half-light diameter of 
the cluster. We used the structural param eters derived from HST 
imaging by Mackev & Gilmore (2003), according to which the 
clusters have half-light radii of r e ff = 8'.'2 (Fornax 3), r e ff = 5'.'5 
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Fig. 2. Signal-to-noise per (re-binned) pixel versus wavelength. 
Below 5200 A, one pixel corresponds to 0.025 A; above 5200 A 
a pixel spans 0.030 A. 

(Fornax 4) and r e ff = 9? 6 (Fornax 5). Two 2400 s exposures 
were obtained in each direction, leading to a total (science) ex- 
posure time of 9600 s for each cluster. FigureQ]shows the acqui- 
sition image for one of the Fornax 4 spectra. Even this cluster, 
located near the centre of the Fornax dSph, clearly stands out 
above the background. However, because of the limited length 
of the UVES slit (10") separate sky exposures were required to 
properly subtract the sky background from the science spectra. 
These were obtained by scanning the slit parallel to the science 
exposures, but at offsets of 2 '5 from each cluster. Each sky scan 
had an exposure time of 1200 s, and two such scans bracketing 
each science scan were co-added and used as a sky reference for 
the science exposures. We verified that these sky exposures were 
made at locations where no bright, individual stars were evident. 

Initial reduction of the spectra was carried out with version 
3.3.1 of the standard UVES pipeline provided by ESO, running 
within the Gasgano environment. The pipeline performs bias- 
and inter-order background subtraction and flat-fields the spec- 
tra. It then carries out wavelength calibration (using calibration 
lamp spectra) and rebins the echelle orders to a common spatial 
scale. Finally, the echelle orders are merged and a single two- 
dimensional, wavelength calibrated spectrum is output for each 
of the two CCD detectors. Using standard tools in the IRAlfl 
package, we then co-added the two sky exposures correspond- 
ing to each science exposure, subtracted the sky from the sci- 
ence data, and then co-added the four individual scans for each 
cluster. Prior to the pipeline reduction, minor detector blemishes 
were repaired with the f ixpix task in IRAF, but the red CCD 
has a group of hot columns that remain visible as an artifact that 
repeats every ~ 50 A in the final spectra. The 2-D spectra were fi- 
nally collapsed to one-dimensional spectra and bad regions were 
masked out (assigned zero weight) for further analysis. Fig. [2] 

1 IRAF is distributed by the National Optical Astronomical 
Observatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc. under contract with the National Science 
Foundation 




Fig. 3. An empirical luminosity fu nction for the globu lar cluster 
M92, compared with a model from lDotter et alj (120071) . 

shows the signal-to-noise ratio (averaged over 200 A intervals) 
as a function of wavelength for each cluster. The discontinuity at 
5200 A corresponds to the transition from the "blue" (EEV) to 
the "red" (MIT-LL) CCD detector. 

3.2. Colour-magnitude diagrams and luminosity functions 

To derive stellar parameters for Fornax 3 and Fo r nax 5 we used 
the HST-based photometry of iBuonanno et al ] (119981) . kindly 
provided to us by C. E. Co rsi (priv. comm . ). For Fornax 4 
we used the photometry from lBuonanno et"al] (U999), available 
through the CDS archival We transformed the photometry in 
the instrumental (WFPC2) system to standard J ohnson-Cousins 
V and / magnitudes using the relations in Holt zman et al.l 
(Il995h . The colour-magnitude diagrams (CMDs) of Fornax 3 
and Fornax 5 reach a limit of My « +5, just below the main se- 
quence turn-off. Still, stars below this limit contribute about 25% 
of the V-band luminosity (for the luminosity functions shown in 
Fig-El see discussion below), a non-negligible fraction. The data 
are somewhat shallower for Fornax 4 (reaching to My x +4) 
which is located in the more crowded central regions of the 
galaxy and was imaged on the WFPC2 WF3 chip (with a pixel 
scale of 07 10), rather than the PC detector (07045). We extrapo- 
lated below the observational limits using theoretical isochrones 
that have been extensively tes ted against deep H ST observations 
of Galactic globular clusters dDotter et al.ll2007l) . 

To properly model the integrated light, we need to know not 
only the shape of the CMDs, but also the relative numbers of 
stars populating each region of the CMDs. Here, it was less 
straight-forward to use the HST photometry since it is based 
on a combination of two datasets (a "deep" and a "shallow" 
exposure) for each cluster that cover different radial ranges. In 
addition, no detailed information is available about the com- 
pleteness of the photometry. For the luminosity functions (LFs) 
we therefore relied more heavily on the theoretical models of 
iDotter et"al] J2007). Fig.[3]shows a comparison of empirical and 

2 http://cdsarc.u-strasbg.fr/ 
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Fig. 4. Binned colour-magnitude diagrams used in the generation of synthetic integrated-light spectra. The plus markers denote 
individual stars, filled circles are bins of mult i ple sta rs, and open squares are model extrapolations from the iDotter etal] (120071) 
isochrones. For Fornax 4, the iBuonanno et all ([l999) ridge line was used near the main seq uence turn-off and for the subgiants, 
shown with filled triangles. Theoretical isochrones used to model the luminosity functions (IDotter et alj|2007l) are shown as red 



curves. 



model LFs for the metal-poor Galactic globular cluster M92. The 
empiri cal L F was obtained by combining the data of iLee et"al1 
d2003l) a nd Paust et all d2007l) . while the model LF assumes a 
Salpeter (Il955h stellar mass function (i.e. a uniform power-law 
dN/dM cc M~ 235 ). Overall, the agreement is very good within 
the range of absolute magnitudes for which empirical data are 
available (i.e., down to My = +9). For stars brighter than the 
main sequence turn-off, the LF shape is mostly determined by 
stellar evolution and does not depend significantly on the as- 
sumed stellar mass function. At masses below ~ 0.5M Q , cor- 
responding to My ^ +8, the mass function in globular clusters 
(and hence the LF) is in general significantly shallower than the 
Salpeter law (Bastian et al J 120101 : iParesce &De MarchJ koOO: 
lDeMarchietal.ll2010h . so an LF based on a pure Salpeter law 
will overestimate the contribution from stars below this limit. 
The detailed shape of the low-mass end of the MF might be 
estimated if the degree of dyn amical evolution experienced by 
each cluster were known (e.g. [Kruiissenll2009l) . but this is be- 
yond the scope of this work. Given the uncer tainty about t he MF 
at the low-mass end, we simply adopted the iDotter et alJ d2007l) 
LF down to a magnitude limit of My — h9 and ignored fainter 
stars. Stars below this limit do contribute to the mass (and hence 
affect the M/L ratio), but their contribution to the integrated light 
is small, less than 5% (in the V-band) if the Salpeter slope ex- 
tends down to the lowest masses, and probably less than that for 
any real cluster. Nevertheless, the behaviour of the low-mass end 
of the LF is a significant contributor to the overall uncertainties 
on our abundance determinations for some elements (Sect. [33]). 

It is not desirable to rely fully on theoretical LFs, particu- 
larly for the horizontal branch (HB). We therefore adopted a hy- 
brid approach in which we based our modelling of the integrated 
light purely on the observed CMDs for the brighter stars (includ- 
ing the HB), combined with the theoretical LFs for fainter stars. 
For Fornax 3 and 5 we used LFs computed for [Fe/H]=-2.2, an 
age of 13 Gyr and [a/Fe]=+0.2, while we assumed [Fe/H]=-1.5, 
an age of 12 Gyr and [a/Fe]=0.0 for Fornax 4. It is important 
to note that the switch from empirical information to theoretical 
modelling occurs at a significantly brighter threshold for the LFs 
than for the overall CMD shape. Specifically, we used empirical 
LF information for red giant branch (RGB) stars brighter than 
My = and for the HB, while the Dotter et al. LFs were used 
elsewhere. These limits are, of course, specific to the dataset 



used here, and would be different in other applications of our 
technique. Note that the stars used in the CMD modelling were 
not selected to fall strictly within the area covered by the slit 
scans. The photometry extends out to larger radii (~ 15") than 
those included in the scans, but does not include stars closer 
than ~ 2" to the centre. Our modelling of the integrated light 
is therefore subject to stochastic (and possibly systematic) varia- 
tions in the number of stars of different types within the slit scan 
areas. The potential impact of these effects is discussed in detail 
in Sect. 



In some globular clusters, blue stragglers (BSs) contribute 
significantly to the integrated light. The CMDs of the three GCs 
studied here do not show any significant population of blue 
stragglers, and they are not included in the modelling. Inclusion 
of BSs would, in principle, be straight forward for clusters where 
the location of the BSs in the H-R diagram and their relative 
numbers are known. However, for clusters with significant BS 
populations a potential comp lication is that B Ss may exhibit pe- 
culiar abundance patterns dLovisi et alj|2012l) . 

The foreground extinction towards the Fornax dSph is small 
but non-negligible. For Fornax 3 and Fornax 5 we assumed 
E(V-I) = 0.05 mag and E(V -I) — 0.08 mag and a d istance mod- 
ulus (m - M)o = 20.68 mag (Buo nanno et alll998l) . For Fornax 
4 ther e is a larger range of E (V-I) values quoted in the liter- 
ature. iBuonanno et alJ dT999) find a relative ly large foreground 
reddening of E(V-I) = 0.15 mag, while iGreco et all (120071) 
find E(B-V) = 0.08 ± 0.03 from the blue edge of the RR Lyrae 
strip, corresponding to E(V-I) = 0.093. For the larger reddening 
value of E(V-I) = 0.15, the RGB ridge line of Fornax 4 becomes 
very simi lar to that of the o ther Fornax GCs. The smaller E(V-I) 
value of Greco et al ] d2007h leads to an intrinsically redder RGB, 
more consistent with a higher metallicity. In this paper we will 
adopt E(V-I) - 0.09 for Fornax 4, along with the same distance 
modulus as for the other two clusters. Varying the distance mod- 
ulus within the range found by various authors (~ 0.1 mag) has a 
negligible effect on our results. We will discuss the uncertainties 
related to the extinction correction in more detail below. 

Tables [T][3] (available in the on-line edition only) list our fi- 
nal adopted CMDs for each cluster. The first two columns list 
the My magnitude and V — I colour for each cmd-bin, while 
the following column lists the weight of the corresponding bin 
when computing the integrated spectra. We used the photometry 
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for the brightest 20 RGB stars directly, and the corresponding 
weight for each bin is then simply one. Elsewhere in the CMDs, 
larger numbers of stars were averaged together per bin, rang- 
ing from 5 stars per bin just below the 20 brightest individual 
stars to 100-150 stars per bin at the faint end. Each CMD then 
consists of about 100 bins. An "E" flag after the weight indi- 
cates that the weight is simply the number of observed stars per 
bin. An "S" flag instead indicates that the weights were deter- 
mined from the model LFs, normalized to the empirical ones 
over the range > My > -2. The 5th column indicates the 
cumulative V-band flux for each bin. It can be seen that the 
HB contributes ~5%-20% of the V-band flux, and a significant 
fraction of the flux comes from stars below the observational 
limit. Finally, the last column indicates the nature of the stars in 
each bin. Here, "IS O" indicates that the colours come from the 
iDotter et all d2007h isochrones, while "RIDGE" (in th e case of 
Fornax 4) is the ridge line from lBuonanno et"al] d 1999b . 

The binned CMDs are shown in Fig. [4] Note the smooth con- 
tinuation from empirical to theoretical CMDs somewhat below 
the main seque nce turn-off. For Forna x 4 we used the empirical 
ridge line from lBuonanno et all ([1999) in the range 2 < My < 4 
where the contamination of the CMD from the general field 
beco mes significant. We have drawn the theoretical isochrones 
from lDotter etalJ d2007l) on top of each CMD. They fit the ob- 
served CMDs fairly well, although some slight mismatches are 
seen at the faint end of the RGB. 

4. Abundance analysis from integrated spectra 

We now explain in more detail our approach to measuring de- 
tailed chemical abundances from the integrated-light spectra. In 
essence, we simply calculate a series of SSP models at high 
spectral resolution, taking into account all stars across the H- 
R diagram and making use of model atmospheres and spectral 
synthesis to compute the individual model spectra for a given 
chemical composition. We calculate these models "on-the-fly" 
while iteratively solving for the best-fitting abundances, i.e., we 
do not make use of pre-computed spectral libraries. Naturally, 
some trade-offs are involved. To make this computationally fea- 
sible, we rely (at least for the time being) on classical model at- 
mospheres and synthetic spectral calculations, sp ecifically those 
computed wit h the publicly available ATLAS9 dKurucz! 1970) 
and SYNTHE dKurucz & Furenlidlll979l: iKurucz & Avretflll981l) 
codes. We use the versions of these cod es available as of Apri l 
2011 from the website of F. CastellQ dSbordone et al.ll2004|) . 
from which we also obtained ancillary data (line lists, molec- 
ular data, opacity distribution functions, etc.). The Kurucz codes 
were compiled with the Intel Fortran compiler for Mac OS 
X. The Kurucz models are based on the common assumption 
of Local Thermodynamic Equilibrium (LTE), as well as being 
one-dimensional, static and plane parallel. Although none of 
these approximations are ever completely valid, they simplify 
the calculations enormously and are often made in the anal- 
ysis of globular cluster stars. In younger clusters, where stars 
with higher temperatures and/or low surface gravities contribute 
a significant fraction of the light, non-LTE effects will likely 
start to play a more significant role. The limitations associated 
with the use of clas sical model atm ospheres have been thor- 
oughly discussed by Asplundj d2005l) . The line lists are based 
on those available from the Kuruc z web sit^l with some mod- 
ifications (Castelli & Hubrig 2004). They include essentially all 

3 http://wwwuser.oat.ts.astro.it/castelli/ 

4 http://kurucz.harvard.edu/linelists.html 



atomic transitions with laboratory data, as well as several di- 
atomic molecules (CH, MgH, NH, OH, SiH, H 2 , C 2 , CN, CO, 
SiO and TiO) and H2O. Of the elements analysed here, hy- 
perfine splitting is included for Sen (/15527,5658 A), Mnn 
(/14756,4765,4783,4824 a) and Ban (i4554, 4934, 5854, 6142 
A). The line list includes damping constants for radiative, Stark 
and van der Waals damping, along with literature references for 
these constants. When literature values are not available, they are 
computed internally in S YNTHE using an approximate expression 
dKurucz & Avrettll 198 ll) . 

4.1. General overview 

In detail, then, we proceeded as follows: We first used the CMDs 
in Fig. |4] to derive physical parameters for each cmd-bin (where 
a "cmd-bin" may refer to a single star or a bin containing several 
stars). From the mean V-I colours of each cmd-bin we derived 
effective temperatures T e s and bolometric corrections BC, using 
the Kurucz colour- T e g transformations. Surface gravities, logg, 
then followed from the standard relation 

log^log, + log[(^)(^)(^) ] (1) 

Here we assumed masses of 0.9 M for the empirical data, cor- 
responding roughly to the main sequence turn-off mass for a 10 
Gyr-old population, while the masses tabulated in the isochrones 
were used for the theoretical data. At the level of precision to 
which the logg values are required, this is sufficiently accu- 
rate. Because the transformations from colour to T e ff and BC 
depend (weakly) on logg, we iterated a few times on Eq. dTJ. 
As the transformations also depend on metallicity, we started 
the analysis with initial "g uesses" based on the RGB slope 
dBuonanno et al.ll998l ll999). but then self-consistently used our 
own metallicity determinations. 

Having established the physical parameters, we com- 
puted an ATLAS9 model atmosphere for each cmd-bin. We 
used the NEWODF opacity distribution functions (ODFs; 
ICastelli & Kuruczl l2003) with solar-scaled abundances, but ver- 
ified that the results did not change significantly for a-enhanced 
ODFs (Sect. 15.31 1. We then computed synthetic spectra with 
SYNTHE. In order to speed up this procedure, the TiO lines were 
skipped for stars w ith > 4500 K wher e the TiO features 
vanish anyway (e.g. lMould & McElrovll 19781) . The spectra were 
initially computed at very high resolution (R = 500000). The 
synthetic spectra were then scaled to the stellar luminosity cor- 
responding to each cmd-bin, multiplied by the weight of the 
cmd-bin, and finally co-added to a single integrated spectrum. 
The integrated model spectra were broadened to match the data, 
and the ratio of the continua of the model- and synthetic spec- 
tra were fitted with a 3rd order spline function or a polynomial, 
depending on the wavelength range fitted. The broadened model 
spectra were then multiplied by the spline/polynomial fit, sub- 
tracted from the observed spectra, and the x 1 evaluated. We in- 
clude an option to assign different weights to different parts of 
the spectra on a pixel-to-pixel basis. In practice, only weights of 
1.0 (good) or 0.0 (bad) were used but intermediate values are in 
principle also allowed. Note that we make no formal distinction 
between "continuum" or "feature" regions for the derivation of 
abundances; the entire spectrum is fitted. Prior to this analysis, 
the observed spectra were shifted to zero radial velocity by com- 
parison with a model spectrum. 

The best-fitting abundances were determined by iteratively 
searching for the abundances that minimized the x 1 - In prac- 
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tice, this was implemented by calling ATLAS9 and SYNTHE from 
a "wrapper" code written in the Python language. This code 
allows us to fit for arbitrary combinations of individual abun- 
dances simultaneously. Elements can be forced to vary in lock- 
step (e.g., the a-elements), while others may be kept fixed at 
user-specified abundances (e.g. from a previous run). The overall 
scaling of all abundances (relative to Solar) can be kept fixed, or 
allowed to vary as a free parameter during the fit. The smoothing 
and microturbulence can also be kept fixed, or allowed to vary 
as free parameters. When multiple parameters a re fit simultane- 
ously we use the "down hill simplex" algorithm dNelder & M ead 
1 965; iPressetaill 19921) . while a "golden section" search is used 
when only one abundance is being fit. In the latter case, the code 
can also output one-sigma error intervals (by varying the abun- 
dance until x 1 - Xmm + 1)' based on the errors provided for the 
observed spectrum. Abundances quoted in this paper are given 
relative to the Solar composition (iGrevesse & Sauvall [T998). 



The model atmospheres were only computed once, at the be- 
ginning of each run. Computing new atmospheres each itera- 
tion would significantly increase the computational burden, and 
we instead leave it up to the user to provide reasonable initial 
guesses for the chemical composition of the stars. If these are 
poorly known, it may be necessary to iterate once or twice. To 
make efficient use of modern multiple-core computers, we have 
set up the Python code to run multiple instances of ATLAS9 
and SYNTHE in parallel. The analysis in this paper was carried 
out on a 12-core Mac Pro workstation. The time required per 
fit depended on many factors, ranging from a few minutes (if 
only a single abundance was measured over a small wavelength 
range, and no errors were required) to a few hours (fitting mul- 
tiple abundances over a 200 A region). We note that while we 
are currently relying on the Kurucz codes as a "backbone", most 
of the Python code is not specific to ATLAS9 and SYNTHE and 
implementing other model atmospheres and/or spectral synthe- 
sis codes should be relatively straight forward. In parti cular, 
the co de currently contains "hooks" for the MOOG code dSnedenl 
Il973h although this is not yet fully implemented. 



4.2. Microturbulence 

The microturbulent velocity v t can be fitted as a free parame- 
ter (which is then, in the current implementation, common to all 
stars) or specified individually for each bin in the input CMD. 
We carried out a set of test runs in which we fitted for the 
overall scaling of the abundances, the Fe and a-element abun- 
dances, and the microturbulence, in adjacent 200 A bins from 
4200 A - 6200 A. For the more metal-poor clusters (Fornax 3 
and Fornax 5) we found v, to be poorly constrained by the data, 
with large bin-to-bin dispersion. This is not surprising, given 
that the equivalent widths of the mostly weak lines in these 
spectra are relatively insensitive to the amount of microturbu- 
lence. For Fornax 4, we instead found the microturbulence to 
be well constrained with a mean (y t ) = 2.0 km s _1 with a bin- 
to-bin dispersion of cr(v,) = 0.18 km s -1 . The use of a single 
v t value for all stars in a cluster is probably an oversimplifica- 
tion, however. While our best-fit value for Fornax 4 is typical of 
the values de termined for individual bright RGB stars in glob- 
ular clusters dKraft et al.lll992t IPilachowski et alll 19961) includ- 
ing the Fornax GCs ( Letarte et al.ll2006l) . v> tends to decrease 
with increasing surface gravity. For subgia nts and stars near th e 
turn-off, a typical value i s v, ~ 1 km s dCarretta et al.l l2004). 
iMcW illiam & Bernstein! d2008l) parameterized this variation as a 
linear function of log g, constrained by the log g and v, values for 
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Arcturus and the Sun. For our purpose, it was more convenient 
to parameterize the trend as a function of direct observables, 
since this allowed us to assign the v< values to each cmd-bin in a 
transparent and straight forward way, prior to further analysis. A 
plot of our derived logg values vs. My showed these two quan- 
tities to be well approximated by a linear relation for My < +4. 
We therefore adopted two reference points: (logg, v,)\ = (1.0, 
2.0 km s _1 ) and (logg, v t )2 = (4-0, 1.0 km s _1 ), corresponding 
to (M v ,v,)i = (-1.81 mag, 2.0 km s" 1 ) and (M v ,v r ) 2 = (+3.64 
mag, 1.0 km s ). For cmd-bins brighter than {My)i that repre- 
sented non-HB stars, we assigned v, values based on a linear fit 
to these two points. For cmd-bins fainter than (My)2 we assumed 
a constant v, = 1 km s , and for cmd-bins rep resenting horizon- 
tal br anch stars we assumed v, = 1.8 km s _1 (IPilachowski et al.l 
119961) . 

4.3. Smoothing of the model spectra 

The model spectra must be degraded to fit the data. This is 
mainly due to the finite instrumental resolution and internal ve- 
locity dispersions in the clusters, although other effects also con- 
tribute to the broadening of spectral features (e.g. stellar rotation, 
macroturbulence, etc.). These effects are all expected to scale as 
0"smoothM - const ((Tsmooth being the dispersion in wavelength 
units of the Gaussian kernel used for the smoothing). This is 
clearly the case for the various velocity -related broadenings, but, 
in our case, also (roughly) for the instrumental broadening which 
is dominated by the slit width: Ignoring distortions by the in- 
strument optics, the projected slit image has a roughly constant 
width in pixel units, s p , set mainly by the focal length of the 
camera. From the grating equation, the dispersion scales with 
echelle order m as dx/dA cc m, so that the width of the slit will 
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also be roughly constant in wavelength units (s^) within a given 
echelle order. However, since the product mA ~ const across the 
different echelle orders, we expect the ratio s^/A - const over a 
wide wavelength range, as well. Not accounting for other effects, 
the instrumental resolution s^/A should thus be roughly constant, 
modulated by a trend with wavelength within each echelle order. 
This latter effect is small, however, because the UVES echelle 
grating is used in very high orders (m = 100 - 145) so that 
the intra-order variation in the ratio sx/A is only expected to be 
~ 1 %. We therefore ignored this more complicated behaviour, es- 
pecially since the instrumental broadening turned out to be rela- 
tively small compared to that introduced by the relative motions 
of the stars in the clusters, as discussed below. 

While our fitting procedure allowed us to solve for <x smoot h as 
a free parameter, this is not generally desirable as it introduces 
an extra degree of freedom in the fits that in some cases may 
dominate the uncertainties on the derived abundances. We there- 
fore carried out a new set of fits to each spectrum in 200 A bins, 
now fixing the microturbulence v, to the relation discussed in 
Sect. 14.21 but still allowing cr smooth and the Fe and a-element 
abundances to vary. These runs then provided the smoothing 
and initial estimates of chemical composition used in subsequent 
analysis. Figure [5] shows the resulting cr smoo th/A vs. A relations 
for each cluster. While cr smoot h/A is roughly constant with A, as 
expected, the residuals in each bin do display significant system- 
atics when comparing the different clusters (some bins are sys- 
tematically high, others systematically low). This serves to em- 
phasize that significant systematic errors might occur if <x smoot h 
were allowed to vary freely. 

From this initial set of fits, we also obtained [Fe/H] = -2.35 
(Fornax 3), [Fe/H] = -1.40 (Fornax 4) and [Fe/H] = -2.06 
(Fornax 5). The bin-to-bin scatter is o-[p e /H] = 0.12 dex, 0.06 
dex and 0.16 dex for the three clusters, respectively. These fits 
also returned mean [ff/Fe] ratios of [a/Fe] = +0.12 (Fornax 3), 
[ar/Fe] = -0.02 (Fornax 4) and [or/Fe] = 0.14 (Fornax 5) with 
bin-to-bin scatter of cr^/pe] =0.11 dex, 0.07 dex and 0.18 dex. 
The smaller bin-to-bin scatter for Fornax 4 is presumably due 
to the stronger features in the spectrum of this cluster. While 
these numbers are subject to further refinement and more care- 
ful analysis in the following sections, we already note that our 
[Fe/H] measurement for Fornax 3 agrees re assuringly well with 
the value of [Fe/H] = -2.4 ± 0. 1 derived bv lLetarte etal] (l2006h 
from high-dispersion spectroscopy of individual stars. It is also 
clear that Fornax 4 is significantly more metal-rich than Fornax 3 
and Fornax 5, and may have a somewhat lower [a/Fe] ratio than 
the other two clusters. 



5. Results and uncertainties 

Our initial fits reported in Sect. [4] were carried out in fixed 
200 A wavelength intervals. However, the measurable features 
of individual elements are not uniformly distributed across the 
spectra, and for our more detailed abundance analysis we con- 
centrated on a number of selected wavelength windows. These 
were defined by visually inspecting the spectra with individual 
absorption features labelled. Although the higher metallicity of 
Fornax 4 in principle allowed the measurement of a larger num- 
ber of weaker lines, we kept the same windows for all three clus- 
ters for consistency reasons. While our procedure allowed us to 
fit for multiple abundances simultaneously, in practice it sufficed 
to fit one element at the time in most cases. 

We first fit for the Fe abundances. Because of the large num- 
ber of Fe lines found throughout the spectra, we simply divided 



the spectra into 200 A bins as before and fit for the Fe abun- 
dance in each of the resulting 10 bins while keeping the smooth- 
ing in each bin fixed according to the relations in Fig. [5] The Fe 
abundances were obtained by fitting not only for Fe, but also for 
other elements that contained prominent absorption lines in the 
respective bins. However, for most bins it made very little dif- 
ference whether we actually let these other abundances vary or 
kept them fixed. The notable exception was the 4200 A - 4400 A 
range which contains prominent CH absorption bands; if the C 
abundance was forced to solar-scaled values this generally led to 
significantly higher [Fe/H] values in this wavelength range (by 
up to 0.4 dex in the most extreme case of Fornax 4) compared to 
the other bins. By letting the C abundance vary, the [Fe/H] abun- 
dance in the 4200 A - 4400 A range was instead consistent with 
those obtained at other wavelengths. From the fits to the 4200 A 
- 4400 A and 4400 A - 4600 A bins, we found C to be signif- 
icantly depleted compared to Fe, with [C/Fe] ratios in the range 
-0.5 to -0.7. It is well known that the [C/Fe] ratio decreases as a 
function of luminosity for RGB stars, particularly at low metal- 
licities whe re C can be depleted by nearly 1 dex for the brightest 
RGB stars (iMartell et al.ll2.Q08b . In our integrated-light spectra 
we only measured the luminosity-weighted mean [C/Fe] values, 
which are clearly significantly sub-solar. However, it should be 
kept in mind that these values were obtained by scaling the C 
abundance of all stars in the H-R diagrams, and as such our in- 
tegrated synthetic spectra should be considered less reliable in 
the 4200 A - 4400 A region where the CH features are strong. It 
should, in principle, be possible to allow for a (suitably parame- 
terized) variation in the [C/Fe] ratio along the RGB. To be fully 
realistic, such modelling would also need to account for vari- 
ations in the carbon isotope ratios as a function of stellar type. 
The Kurucz code uses a hard-coded ratio of 12 C/ 13 C = 89, i.e. the 
Solar value (e.g. lHarris et al.|[l987h . while giant s tars generally 
have much lower isotop ic ratios of 12 C/ 13 C ^10 dGratton et al.l 
2000; [Keller et al.l l2~001). The integrated-light measurements of 
C likely reflect a combination of effects related to both stel- 
lar evolution and the overall chemistry of the clusters, and are 
clearly complicated to interpret. For these reasons, we do not in- 
clude C in our further discussion of individual elements below. 

When measuring the remaining elements, we kept the overall 
scalings of the abundances fixed according to the derived [Fe/H] 
and [C/Fe] abundances. We then proceeded by first fitting the el- 
ements with more absorption features, starting with Ti and then 
Ca. Abundances of these elements were then fixed at the fitted 
values, and further elements were fitted, etc. Example fits are 
shown in Fig. [6] where we plot the data together with the best- 
fitting model spectra and models where the Fe abundances were 
changed by ±0.3 dex. Both models and data have been smoothed 
with a 5 pixels wide boxcar filter for the purpose of illustration. 
The difference between the spectra of the very metal-poor clus- 
ters Fornax 3 and Fornax 5 and the significantly more metal-rich 
Fornax 4 is clearly seen. 

Our individual abundance measurements are reported in 
Table|4]together with the associated formal errors. When a wave- 
length range is given in the first column, this indicates that sev- 
eral lines are present within that range but the overall abundance 
was obtained from one fit. When instead a single wavelength is 
given, only one absorption line from the element is present, and 
the fit was carried out over a range of +5 A (roughly) centered 
on the wavelength of the feature. For windows wider than 10 A 
the continua of the model- and observed spectra were matched 
using a 3rd order spline with typically 5 knots, while a 2nd order 
polynomial was used for the smaller wavelength bins. 
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Table 4. Individual abundance measurements. 
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.031 


+0 


,346: 


= 0. 


066 


4352 


+0 


.137 = 


= 


.086 


-0 


.344 = 


= 


.086 


+0.526 ± 


0.131 


4353-4400 


-0 


.214 


±0 


,070 


-0 


.364: 


= 


,051 


-0 


.144: 


= 0. 


096 


4571 


+0 


.076: 


= 


.086 


-0 


.044 = 


= 


.066 


+0.307 ± 


0.115 


4520-4660 


-0 


.205 


±0 


,035 


-0 


.073: 


= 


,016 


-0 


.044 d 


= 0. 


041 


4703 


-0 


.195 = 


= 


.056 


+0 


.015: 


= 


.041 


-0.034 + 


0.081 


5235-5330 


-0 


,244 


±0 


,061 


-0 


.034: 


= 


,025 


-0 


.093: 


= 0. 


060 


5167 


-0 


.195 = 


= 


.076 


-0 


.115: 


= 


.051 


+0.137 + 


0.111 


5342-5351 


+0 


.026 


±0 


.076 


+0 


,127: 


= 


,051 


+0 


.307: 


= 0. 


101 


5528 


-0 


.034 = 


= 


.061 


+0 


.046 = 


= 


.056 


-0.034 ± 


0.096 


5410 


+0 


.005 


±0 


.086 


+0 


.076: 


= 


,070 


-0 


,324: 


= 0. 


146 


5711 


+0 


.395 = 


= 


.110 


-0 


.024 = 


= 


.051 


+0.297 ± 


0.151 


[Mn/Fe] 


























[Ca/Fe] 






















4750-4790 


-0 


.523 


±0 


.070 


-0 


,263: 


= 


.025 


-0 


413 d 


= 0. 


.070 


4227 


+0 


.166 = 


= 


.040 


-0 


.003: 


= 


.031 


+0.256 + 


0.061 


[Cu/Fe] 


























4280-4320 


+0 


.147 = 


= 


.036 


+0 


.236 = 


= 


.041 


+0.376 ± 


0.056 


5106 




< 


-1 




-0 


.834: 


= 


.076 


-0 


.555 : 


= 0. 


.210 


4420-4460 


+0 


.266 = 


= 


.025 


+0 


.116 = 


= 


.025 


+0.147 ± 


0.051 


[Y/Fe] 


























4575^1591 


-0 


.134 = 


= 


.096 


-0 


.014: 


= 


.041 


-0.244 + 


0.131 


4355-4425 


-0 


,173 


±0 


.071 


-0 


.289: 


= 


,065 


-0 


,310: 


= 0. 


121 


4878 


+0 


.616 = 


= 


.076 


+0 


.567: 


= 


.070 


+0.607 ± 


0.091 


4884 


-0 


.115 


±0 


.116 


-0 


.334: 


= 


.076 


-0 


.164 d 


= 0. 


171 


5259-5268 


+0 


.366 = 


= 


.066 


+0 


.256: 


= 


.041 


+0.395 ± 


0.086 


5088 


+0 


.266 


±0 


.086 


-0 


.134: 


= 


.070 


-0 


,243: 


= 0. 


146 


5580-5610 


+0 


.116 = 


= 


.041 


+0 


.207 = 


= 


.026 


+0.307 ± 


0.046 


[Ba/Fe] 


























U 1 \J\J — U 1 / J 




.326 = 


= 


.025 




.105: 


= 


.016 


+0.236 d 


= 0.03 


4SS4 




,856 


±0 


.035 




.366: 


= 


,030 




.166: 


= 0. 


.080 


[Sc/Fe] 






















4934 


+0 


,997 


±0 


.056 


+0 


.806: 


= 


,041 


+0 


485: 


= 0. 


.075 


4290-4330 


+0 


.046 = 


= 


.046 


+0 


.127 = 


= 


.046 


+0.026 ± 


0.076 


5846 


+0 


,376 


±0 


.086 


+0 


.185: 


= 


,081 


-0 


,173: 


= 0. 


150 


4350-4440 


+0 


.297 = 


= 


.036 


-0 


.034 = 


= 


.036 


+0.046 ± 


0.061 


6142 


+0 


.707 


±0 


.056 


+0 


.227: 


= 


,051 


-0 


.205: 


= 0. 


121 


4670 


-0 


.454 = 


= 


.156 


-0 


.364 = 


= 


.086 


-0.383 ± 


0.156 


[Eu/Fe] 


























5031 


-0 


.183 = 


= 


.091 


-0 


.034: 


= 


.046 


-0.723 + 


0.150 


4205 


+1 


.527 


±0 


,131 


+2 


,197: 


= 


,096 


+1 


.286: 


= 0. 


.206 


5527 


+0 


.276 = 


= 


.086 


-0 


.124: 


= 


.070 


-0.183 ± 


0.166 


4436 


+1 


477 


±0 


,126 


+0 


.606: 


= 


,146 


+0 


,177: 


= 0. 


.405 


5638-5690 


+0 


.256 = 


= 


.051 


-0 


.064 = 


= 


.030 


+0.105 ± 


0.081 





























Notes. A wavelength range in the first column indicates that this range contains several absorption features from the corresponding element. When 
a single wavelength is given, a single line is fitted within a 10 A window. 



Table 5. Average integrated-light abundance measurements. 





N 


Fornax 3 




Fornax 4 




Fornax 5 








weighted avg 


r.m.s. 


weighted avg 


r.m.s. 


weighted avg 


r.m.s. 


cr vd [km s 1 ] 


10 


7.5 


0.8 


4.5 


0.2 


5.4 


0.6 


[Fe/H] 


10 


-2.33 ± 0.01 


0.09 


-1.42 ±0.01 


0.06 


-2.09 ± 0.01 


0.14 


[Mg/Fe] 


6 


-0.04 + 0.03 


0.22 


-0.04 + 0.02 


0.13 


+0.13 ±0.04 


0.21 


[Ca/Fe] 


8 


+0.25 + 0.01 


0.20 


+0.13 ±0.01 


0.18 


+0.27 ± 0.02 


0.23 


[Sc/Fe] 


6 


+0.18 ±0.02 


0.31 


-0.04 ± 0.02 


0.15 


-0.03 ± 0.04 


0.33 


[Ti/Fe] 


9 


+0.29 ± 0.01 


0.11 


+0.12 ±0.01 


0.15 


+0.24 ± 0.02 


0.13 


[Cr/Fe] 


6 


-0.17 ±0.02 


0.11 


-0.06 + 0.01 


0.16 


+0.02 + 0.03 


0.24 


[Mn/Fe] 


1 


-0.52 + 0.07 




-0.26 ± 0.03 




-0.41 + 0.07 




[Cu/Fe] 


1 


< -1.0(2cr) 




-0.83 ± 0.08 




-0.56 ±0.21 




[Y/Fe] 


3 


-0.02 ± 0.04 


0.19 


-0.25 ± 0.04 


0.09 


-0.26 ± 0.08 


0.06 


[Ba/Fe] 


4 


+0.81 ±0.02 


0.24 


+0.45 ± 0.02 


0.25 


+0.20 + 0.05 


0.31 


[Eu/Fe] 


2 


+ 1.50 ±0.09 


0.03 


+ 1.72 ±0.08 


0.86 


+ 1.06 ±0.18 


0.64 



Notes. A' is the number of individual fits for each element (see Table [4]l. For each entry we give the weighted average of the measurements in 
Tableland the r.m.s. scatter. 



In Table [5] we list the weighted average abundances per el- 
ement for each cluster. In addition to the formal errors on the 
weighted means, we also list the r.m.s. bin-to-bin scatter of the 
individual abundance measurements when applicable. The lat- 
ter may be viewed as a more realistic indicator of the actual 
measurement uncertainties, and may be turned into errors on 
the mean abundances by taking into account the number of in- 



dividual measurements (N), i.e., cr - r.m.s/ VA^ - 1. In most 
cases, this yields errors smaller than 0.1 dex, but nevertheless 
somewhat larger than the formal errors based on the^ 2 analysis. 
That the purely random errors do not fully account for the bin- 
to-bin scatter on the derived abundances can also be seen from 
the fact that some wavelength bins yield systematic higher or 
lower abundances than others. Indeed, the abundance measure- 
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Table 6. Systematic uncertainties on abundance measurements. 





AE(V-1) = 


0.05 


Av, 


= 0.5 km s _J 


M v ,h 


ra 9^13 mag" 




F3 


F4 


F5 


F3 


F4 


F5 


F3 


F4 


F5 


A[Fe/H] 


+0.09 


+0.09 


+0.13 


-0.11 


-0.11 


-0.12 


-0.07 


-0.01 


-0.04 


A[Mg/Fe] 


-0.04 


0.00 


-0.04 


+0.06 


+0.01 


+0.03 


-0.03 


-0.05 


-0.04 


A[Ca/Fe] 


0.00 


+0.01 


-0.02 


+0.03 


-0.04 


+0.03 


-0.10 


-0.10 


-0.10 


A[Sc/Fe] 


-0.02 


-0.05 


-0.01 


+0.01 


-0.07 


+0.00 


+0.09 


+0.03 


+0.05 


A[Ti/Fe] 


-0.02 


-0.02 


+0.03 


-0.01 


-0.06 


-0.01 


+0.00 


-0.04 


-0.02 


A[Cr/Fe] 


+0.02 


0.00 


0.00 


+0.02 


-0.02 


+0.02 


-0.11 


-0.05 


-0.06 


A[Mn/Fe] 


0.00 


0.00 


-0.01 


+0.08 


-0.01 


+0.07 


-0.01 


-0.01 


-0.02 


A[Cu/Fe] 




+0.02 


+0.01 




+0.06 


+0.10 




+0.02 


+0.06 


A[Y/Fe] 


0.00 


-0.06 


0.00 


+0.02 


-0.10 


+0.01 


+0.17 


+0.06 


+0.09 


A[Ba/Fe] 


+0.02 


-0.01 


+0.01 


-0.11 


-0.13 


-0.11 


-0.00 


-0.02 


-0.00 


A[Eu/Fe] 


-0.04 


-0.04 


-0.03 


-0.01 


-0.23 


+0.01 


+0.02 


-0.05 


-0.05 



Notes. (o) Indicates a change in the limiting magnitude of the CMD modelling from M v = +9 mag to My = +13 mag. 



ments depend on other parameters, such as the atomic oscillator 
strengths of which some still remain quite uncertain (see e.g. 
discussion of Ti and Ba in Sect. 15. II ) and hyperfine splitting. 

5.1. Comparison with Letarte et al. results for Fornax 3 

iLetarte et alJ {2006) measured the abundances of several ele- 
ments in individual stars in Fornax 1, Fornax 2 and Fornax 3. 
In the following, we compare our integrated-light results with 
their data for three RGB stars in Fornax 3. 

Fe: L06 find an iron abundance of [Fe/H] = -2.4 + 0.1. This 
is in excellent agreement with our integrated-light measurement 
of [Fe/H] = -2.33. The r.m.s. scatter of 0.08 dex on our [Fe/H] 
measurements corresponds to an error on the mean of ±0.03 dex. 

a-elements: For Ca, our measurement of [Ca/Fe] = +0.25 
agrees very well with the L06 result of [Ca/Fe] = +0.22 + 0.02. 
For Ti, L06 measure [Til/Fe] = -0.06 + 0.05 and [Tin/Fe] = 
0.10 + 0.02. These values are both significantly lower than our 
measurement of [Ti/Fe] = +0.29, which furthermore has a rela- 
tively small r.m.s. scatter of 0.1 1 dex with 9 individual measure- 
ments. As noted by L06, the log gf values for Ti are somewhat 
uncertain, translating into corresponding uncertainties on the de- 
rived abundances. Comparing the log gf values for Ti used by 
L06 with those in the Kurucz line list, we find that the log gf 
values in the Kurucz list are generally lower, with an average 
difference of about A(Kurucz-L06) = -0.15 dex but with a range 
from -0.8 dex to +0.15 dex. Applying the average difference to 
the Ti abundances would bring our measurements in much better 
agreement with those of L06. However, this comparison is only 
indicative as the overlap between the lines we actually measure 
is limited. L06 used a different UVES set-up, covering a redder 
wavelength range that does not include the range from 4200 A - 
4800 A where many of our Ti measurements are made. 

Our [Mg/Fe] ratios are quite low compared to the other a- 
elements, [Ca/Fe] and [Ti/Fe], for all three clusters. Here, the 
comparison with L06 is not very conclusive as L06 find a sub- 
stantial spread in the [Mg/Fe] values for their three stars of 
+0.19 ± 0.09, +0.37 ± 0.07 and -0.35 + 0.08. This spread may be 
related to the antic orrelatio n of Al and Mg that is known to exis t 
in some Milky Way GCs (IShetronelll9"96l iGratton et al.ll200lh . 
Neither we nor L06 measure Al, but it is interesting to note that 
low (solar-like) [Mg/Fe] abundance ratios have al so been found 
from integrated-light spectroscopy of M3 1 GCs dColucci et alJ 



2009). It is therefore clear that caution should be exercised when 
using Mg as an indicator of the [a/Fe] ratio in the integrated light 
of GCs. 

Other elements: Other elements in common between our study 
and that of L06 are the Fe-peak elements Mn, Cr and the heavy 
elements Ba, Eu and Y. For Fornax 3 we find [Cr/Fe] = -0.17 
with an error of ~ 0.05 dex based on the r.m.s. scatter. For their 
three stars in Fornax 3, L06 instead find [Cr/Fe] = -0.28 ± 0.20, 
-0.30 + 0.20 and -0.50 + 0.15. The mean of the three L06 mea- 
surements is -0.39 + 0.10, formally a ~ 2cr difference with re- 
spect to our data. We note, however, that inclusion of low-mass 
stars below the adopted magnitude limit of our CMDs would 
tend to (slightly) decrease the Cr abundances (Sect. 15.31 and 
Table|6]). Furthermore, L06 estimate a 0.21 dex model-dependent 
systematic uncertainty on their Cr abundance. We therefore con- 
clude that our Cr abundance for Fornax 3 is consistent with the 
measurements of L06 within the uncertainties. 

For Mn, we find a significantly sub-solar abundance ratio 
of [Mn/Fe] = -0.52 while L06 find <[Mn/Fe]> = +0.01, in 
both cases with formal error bars less than 0. 1 dex. However, our 
[Mn/Fe] measurement is based on only a single wavelength bin 
in the blue, making it difficult to assess the real uncertainty. This 
wavelength range was not included by L06, while the 6014 A 
and 6022 A lines measured by L06 are too weak in our data to 
provide a reliable measurement. Indeed L06 only quote [Mn/Fe] 
for two or their three stars in Fornax 3, for one of three stars 
in Fornax 2 and for no stars in Fornax 1, attesting to the diffi- 
culty of these measurements. Our lower [Mn/Fe] abundance ra- 
tio is more similar to those seen in metal-poor Milky Way field 
stars and globular clusters and for field stars in the Fornax dSph, 
although [Mn/Fe] ratios closer to Solar have also been derived 
from integrate d-light measurements of LMC globular clusters 
dColucci et al.ll201 2|). We return to this point below. 

Our abundances of Y, Ba and Eu ([Y,Ba,Eu/Fe] = 
-0.02, +0.81, +1.50) are somewhat higher compared to L06 
([Y,Ba,Eu/Fe] = -0.22, +0.22 and +0.90). For Ba in particu- 
lar, the 0.6 dex difference is too large and too significant to be 
explained by measurement uncertainties. Our most discrepant 
measurement (relative to L06) comes from the 4934 A line. This 
line yields systematically higher Ba abundances than other lines 
not only for Fornax 3, but also for the two other clusters in our 
sample. The log gf value for this line appears to be particularly 
uncertain: in the Kurucz line list it h as log g f - -0.458 . L06 
instead use log gf = -0.703, while IColucci et all ((2009) use 
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Fig. 6. Fits to a small part of the UVES spectra (dotted, thick black curves). The thin, solid curves (blue in the on-line edition) show 
the synthetic spectra for the overall best-fitting Fe abundances, while the thin dashed curves (red in the on-line edition) show the 
effect of of varying the Fe abundance by +0.3 dex. The spectra have been smoothed with a 5 pixels wide boxcar filter. 



log gf - -0. 157. For the other Ba lines there are only minor dif- 
ferences between different sources. If we use the higher log gf 
value of IColucci et"ai1 (120091) for the 4934 A line, our Ba abun- 
dance inferred from this line decreases accordingly and is then 
in better agreement with that inferred from the other Ba lines, 
and the average Ba abundances in Table [5] decrease by ~ 0.15 
dex. If we instead use the smaller log gf value from L06, our 
Ba abundance increases correspondingly and the difference with 
respect to L06 becomes even greater. The reason for this dif- 
ference therefore remains unclear. However, some Galactic GCs 
are known to have a substantial internal spread in the heavy re- 
capture elem ent abundances. In M15, the [Ba/Fe] ra tio varies by 
0.5-0.6 dex (ISneden et al.lll997HSobeck et al.ll201 ll) . and is cor- 
related with similar variations in the [Eu/Fe] ratio. These varia- 
tions are seen also in the small sample of Ml 5 stars measure d 
by L06, and in the globular cluster M22 (Mar ino et all 12009). 
The uncertainties on our [Eu/Fe] abundances are large as the 
two lines we can measure are relatively weak and blended with 
several other lines, but it is worth noting that our [Eu/Fe] mea- 
surement does show a similar offset with respect to the L06 data 
as does [Ba/Fe]. This might suggest that the difference between 
our measurements and those of L06 are due to Fornax 3 having 



internally inhomogeneous n-capture element abundances, with 
the three stars measured by L06 happening to have less than 
average abundances of these elements. The light n-capture el- 
ement Y does not ap pear to follow the same trends as Eu and Ba 
(ISobeck et alJl201 1), and the difference between our [Y/Fe] ratio 
and that of L06 is indeed much smaller and not very significant 
given the ~ 0.10 dex uncertainties on both measurements. This 
scenario should, in principle, be verifiable by measuring a larger 
sample of individual stars in Fornax 3. 

5.2. Comparison with the spectrum of Arcturus 

As an additional consistency check, we have applied o ur analysis 
to the spectrum of Arcturus dHinkle & \Vallacel2005l) . The phys- 
ical parameters of Arcturus are well established, and as a starting 
po int for our analysis we used the s tellar parameters determined 
bv lRarmrez & Allende Pr"ietol(l201 ll hereafter Rl 1): T eff = 4286 
K and logg = 1.66. These authors also find an iron abundance 
of [Fe/H] = -0.52 + 0.04, microturbulent velocity v, = 1.74 km 
s and an absolute magnitude My = -0.313 + 0.016. For this 
My, our formula yields v, = 1.725 km s _1 , in close agreement 
with the value of Rll. We consider the step of deriving stellar 
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Table 7. Abundances for Arcturus. 





Rll 


This work 


r.m.s. 


[Fe/H] 


-0.52 ± 0.04 


-0.64 


0.04 


[Mg/Fe] 


+0.37 ± 0.03 


+0.22 


0.06 


[Ca/Fe] 


+0.11 ±0.04 


+0.26" 


0.17 


[Sc/Fe] 


+0.21+0.04* 


+0.11 


0.21 


[Ti/Fe] 


+0.23 ± 0.03 c 


+0.37 


0.31 


[Cr/Fe] 


-0.05 ± 0.04 


-0.04 


0.17 



Notes. <D) Rejecting /14878A {b) Average of Sc i and Sc n. {c) Average of 
Ti i and Ti n. 

parameters from integrated photometry a separate problem; this 
will be discussed below. 

No special modifications to our procedure were required in 
order to model a single star - we simply modelled the spec- 
trum of Arcturus as a "cluster" of one star. However, because 
of the higher metallicity of Arcturus compared to the Fornax 
GCs, a much larger fraction of the spectrum is dominated by rel- 
atively strong absorption lines. In a proper analysis, this should 
be taken into account and strongly saturated lines should be 
avoided, but for consistency with the Fornax GC analysis we 
analysed the Arcturus spectrum using the exact same wavelength 
bins. However, due to the stronger lines in the Arcturus spec- 
trum, we found that a simple spline fit to the model/data ra- 
tio, as we did for the Fornax GCs, was significantly affected 
by several strong features in the observed spectrum which had 
no counterparts in the model. As a result, the continuum levels 
in the model spectra were noticeably lowered compared to the 
data, and the code would converge towards a mean iron abun- 
dance of [Fe/H]=s -0.70, about 0.2 dex lower than the liter- 
ature value. We therefore adopted a slight modification to the 
scaling procedure, in which we first identified "continuum" re- 
gions in the observed spectrum where the flux fx > 0.98/ max , 
where / max was determined as a "running maximum flux" in 
5 A bins. We then iteratively rejected pixels where the differ- 
ence (/model -/data) exceeded 2x the standard deviation evaluated 
over all remaining pixels. With this modification to the scaling, 
we found that the continuum levels of the model- and observed 
spectra agreed well, and the average iron abundance increased 
to [Fe/H] =-0.64 with an r.m.s. scatter of 0.04 dex. This is still 
about 0.1 dex lower than the Rl 1 measurement, although lower 
values have been quoted in the literature (e.g. [Fe/H] =-0.68; 
iGriffin & Lvnas-Gravfl999l) . 

Apart from its overall lower metallicity, the abundance pat- 
terns in Arcturus differ from those of the Su n by being enhanced 
in the a-elements (e.g. Peter son et all l993). Abundances for el- 
ements in common between our work and Rll are compared 
in Table [7] Our measurements agree with those of Rll within 
0.10-0.15 dex. From the 4200A - 4400A and 4400 A - 4600 A 

bins we also derive a roughly Solar [C/Fe] ratio ([C/Fe] 0.05). 

This is significantly lower than the Rl 1 measurement of [C/Fe] = 
+0.43±0.07, while lPeterson et all d 19931) quote [C/Fe]=0.0, sim- 
ilar to our measurement. The difference might be related to the 
use of indi vidual C i lines (Rl 1 ) vs. the CH molecular bands (this 
study and Pet erson et alj|1993l) . 

The transformation from broad-band photometry to stellar 
parameters is an important step in our procedure. Systematic 
errors in these transformations would translate to correspond- 
ing systematics in the abundance measurements. While we have 
used Kurucz transformations for solar-scaled abundances, using 
colour transformations based on ff-enhanced model atmospheres 
only changes T e ff by 10-15 K. This is in agreement with the find- 



ings of Cassis i et al ] (12004 . who found that the V — I colour is 
relatively insensitive to a-enhancement (although other colours, 
like B - V and U - B, are more strongly affected). However, 
colour- r e ff transformations differ significantly from one author 
to another. For example, a compa rison of the Kurucz V — I vs. 
r e if transformations with those of Wor thev & Led d201 lb shows 
differences of 50-70 K, with the Kurucz transformations tend- 
ing to give cooler T e s for giants and hotter r e ff for subgiants. 
Assuming T e s = 4350 K for Arcturus (i.e., an increase of 64 K) 
and repeating the analysis, the iron abundance increases by 0.05 
dex to [Fe/H] = -0.59. Overall, systematic errors in the colour- 
T e s transformations have an effect that is quite similar to that of 
changing the correction for foreground extinction, which will be 
discussed in more detail below. 

5.3. Uncertainties on the abundance measurements 

In Sect.[5]we quantified the uncertainties on our abundance mea- 
surements due to formal errors on individual wavelength bins 
and from the bin-to-bin scatter. However, in this type of analy- 
sis one must also consider a variety of systematic uncertainties, 
some of which we discuss in the following. 

Extinction. An increase in the extinction correction E(V - 1) 
leads to intrinsically bluer colours and higher r e ff for the stars. 
Overall, this has the effect of making most spectral lines in the 
model spectra weaker for a given metallicity, so the abundances 
inferred from the data increase to compensate. In Table|6]we list 
the changes in the abundances for an increase of AE{V-I) = 
+0.05 mag. This corresponds roughly to the range of extinc- 
tion values quoted in the literature for Fornax 4, while it prob- 
ably exceeds the uncertainties for Fornax 3 and 5 by a substan- 
tial margin. For this AE(V -I), the Fe abundances increase by 
A[Fe/H] ~ 0.10 dex. Abundance ratios are much less strongly 
affected, and most of the ratios we measure are hardly affected 
by uncertainties in the extinction correction at this level. The 
most strongly affected ratio is [Mg/Fe], which decreases by 
about 0.03 dex. Increasing the [Mg/Fe] ratios to bring them in 
better agreement with the other a-element ratios would thus re- 
quire a decreased E(V-I), but even if we assume zero extinction 
towards each cluster, this would not make the [Mg/Fe] ratios 
similar to [Ca/Fe] or [Ti/Fe]. 

Microturbulence. An increase in the microturbulent velocity v t 
has the effect of de-saturating strong spectral lines, so a given 
observed line strength corresponds to a lower abundance. The 
magnitude of this effect depends on the strengths of the specific 
lines, with stronger (more saturated) lines being more heavily 
affected. The sensitivity of our results to changes in v, was as- 
sessed by adding an offset of Av, = 0.5 km s _1 to the assumed 
v, values and repeating the analysis. As can be seen in Table [6] 
the net effect of this (rather large) change is to decrease the Fe 
abundances by about 0. 1 dex. Again, abundance ratios are gen- 
erally less affected. The most strongly affected abundance ratio 
is [Ba/Fe]; however, even large changes in v f would not be suf- 
ficient to bring our Ba abundance for Fornax 3 into agreement 
with that of L06. 

Luminosity functions. The LFs are critical ingredients in our 
modelling of the integrated light as they determine the rela- 
tive contributions from stars of different luminosities. Although 
Fig. [3] shows good agreement between the model LFs and em- 
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pirical LFs for M92 within the range of My values for which 
data are available, subtle differences might affect the abundance 
determinations at some level. To quantify this, we first replaced 
the Dotter et al. theoretical LFs with the empirical M92 LF for 
Fornax 3 and 5. For Fornax 4 we used an emp irical LF based 
on the globular cluster M12 (lHargis et al.ll2004f) . We found that 
[Fe/H] values as well as individual abundance ratios changed by 
at most 0.010-0.015 dex. It thus appears that uncertainties due 
to the LF shape within this magnitude range are relatively minor. 

At fainter magnitudes, dynamical effects become increas- 
ingly important. Our exclusion of stars fainter than My = +9 
effectively corresponds to the extreme case in which no stars 
with M $ QAM Q are present. To gauge the effect of fainter stars, 
we repeated the analysis for a limiting magnitude of My = +13, 
assuming a Salpeter IMF extending down to O.1M . The last 
column in Table [6] shows the effect of adopting this fainter mag- 
nitude limit. The abundances are indeed sensitive to the low- 
luminosity stars, with some abundance ratios being at least as 
sensitive as the overall [Fe/H] abundance. For realistic clusters, 
the assumption of a Salpeter IMF down to 0.1M o will probably 
overestimate the contribution from the fainter stars, so we expect 
the changes to be smaller than the entries in the last columns of 
Table |6l probably on the order of a few times 0.01 mag for most 
ratios. Nevertheless, it is clearly important for high-precision 
abundance analysis to pay attention to even the faintest stars, 
and the high mass-to-light ratios for Fornax 3 and 5 might indi- 
cate the presence of a relatively large number of low-mass stars 
(iDubafh et alJll992l and Sect. 1331). With appropriate data, the 
contribution of dwarf stars to the integrated light might even be 
constrained as part of the analysis its elf via features that are 
strongly sensitive to surface gravity (Mie ske & KroupaH 2008; 
Conrov & van Dokkurrj20li . 



Opacity distribution functions. When computing the ATLAS9 
model atmospheres we made use of pre-computed ODFs with 
solar-scaled abundance patterns. The ODFs are also available for 
c-enhanced abundances ([a/Fe] = +0.4). Computing the model 
atmospheres with the a-enhanced ODFs and redoing the anal- 
ysis, we found that this had a negligible effect on the results: 
[Fe/H] values changed by less than 0.005 dex, and most abun- 
dance ratios by less than 0.01 dex. The most extreme change 
was for the [Eu/Fe] ratio of Fornax 4, which decreased by 0.014 
dex when using a-enhanced ODFs. 



Stochastic sampling of the stellar IMF. It is well known that 
the finite number of stars in a cluster can cause its integrated 
colours to deviate strongly from those pred icted for a continu- 
ously populated stellar mass function (e.g . iGirardi et alJI 1995: 
Bruzuall 12002 : Fouesneau & Lanccinl 12010: Porjescu & Hanson 
20 lOtlSilva- Villa & Larsenll201 ll) . We assessed the effect of ran- 



dom IMF sampling on our abundance determinations by boot- 
strapping. The original CMDs were resampled randomly with 
replacement and all steps of the analysis were then carried out 
again. This procedure was repeated 10 times, and the uncertain- 
ties due to stochastic IMF sampling were then estimated as the 
standard deviation of the 10 resulting abundance determinations. 

Due to the rather time consuming nature of this experiment 
we did not assess the effect on the full range of abundance ra- 
tios, but restricted the analysis to a few cases: [Fe/H], [Ti/Fe], 
[Ca/Fe], [Mg/Fe] and [Ba/Fe]. For [Fe/H] we found uncertain- 
ties of cr[ Fe / H ],stoc = 0.02 dex, 0.03 dex and 0.05 dex for Fornax 3, 
Fornax 4 and Fornax 5, respectively. The scatter on the elemental 
abundance ratios was always less than on the [Fe/H] values. 



Clearly, this procedure only gives an approximate descrip- 
tion of the actual uncertainties. On the one hand, the stars in the 
HST CMDs are actual stars present in each cluster. Therefore, if 
our CMDs corresponded exactly to the rectangular regions cov- 
ered by the drift-scanned spectra, stochastic fluctuations would 
play a less significant role than suggested by the bootstrapping 
experiment. However, even in this case it should be remembered 
that our CMD modelling is fully empirical only for the brighter 
stars. For fainter stars we rely on model LFs to define the weight 
of each cmd-bin. On the other hand, the CMDs do not corre- 
spond exactly to the parts of the clusters for which we have 
spectra. The central regions are too crowded to obtain reliable 
photometry, and there is no photom etry within ~ 2" of the clus - 
ter centres. From the profile fits of Mackev & Gilmore (2003), 
about 10-15% of the flux falls within the central region excluded 
from the photometry. Additionally, the observed CMDs include 
stars from regions of the clusters not covered by the spectra. 
Roughly 2/3 of the stars in the CMDs fall within the slit area, 
while roughly 1/3 fall outside. Given the significant overlap, the 
true impact of stochasticity is likely smaller than estimated by 
our bootstrapping procedure. 



Radial gradients. To test whether radial gradients in the inte- 
grated properties of the clusters contribute significantly to the 
uncertainties, we divided the spectral scans into an inner section 
(the central +4 pixels or ±1'.'8) and an outer section (the outer- 
most 7 pixels on either side of the centre). The inner sections 
approximately span the radial range for which CMD data are 
missing, although the scans of course extend out to the half-light 
radius in both cases. We then analyzed each section separately. 
The differences were found to be very minor. For [Fe/H] we 
found differences (outer-inner) of 0.001 ± 0.009 dex (Fornax 3), 
-0.029 + 0.006 dex (Fornax 4) and -0.001 ± 0.010 dex (Fornax 
5). Considering these very small differences, we did not look at 
other abundance ratios, but it appears safe to conclude that no 
strong radial trends are present. 



Background subtraction The subtraction of background light 
from the spectral scans may also be affected by uncertainties 
due to stochastic sampling of the underlying field population. 
However, as can be seen from Fig. Q] the contrast between the 
GCs and background is quite high even for Fornax 4. In fact, 
with a central surfac e brightness of jxy - 23.3 m ag arcsec~ 2 
for the Fornax dSph dlrwin & Ha tzidimitriou 1995), the back- 
ground is everywhere dominated by the night sky rather than by 
the Fornax dSph itself. The total background signal, including 
night sky, constitutes ~3%-6% of the counts in the science ex- 
posures, so the underlying field population in Fornax contributes 
only marginally to the signal in each spectrum . This is corrobo- 
rated by the analysis of iBuonanno et al.l (1 19991) . who find a ratio 
of cluster vs. field stars of about 4 at about 5" from the centre 
of Fornax 4, correspo nding to the end points of th e slit in Fig.Q] 
Integrating the fits of iMackey & Gilmoji (120031) out to r — 5 
then gives a number ratio of cluster- to field stars of about 9. 
Converting this to a luminosity ratio is less straight forward as 
the field- and cluster CMDs are quite different, due to the ex- 
tended star formation history of the field. However, if we inte- 
grate the central surface brightness of the Fornax dSph over the 
central 5" 5 of Fornax 4, and compare th is with the appar ent inte- 
grated cluster magnitude of V = 13.57 (Webbink 1985) (adding 
0.75 mag to find the magnitude within r e «) we find that the field 
contributes about 2.5% to the signal. We therefore conclude that 
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Table 8. Velocity dispersions, dynamical masses and mass-to- 
light ratios. 





era, (km s 1 ) 


n, (pc) 


M vir (10 3 M o ) 


My 


M/L v 


For 3 


6.5 ± 0.2 


7.2 


5.2 + 0.4 


-8.2 


3.5 ±0.3 


For 4 


4.1 ±0.1 


4.9 


1.4 ±0.1 


-7.4 


2.0 ±0.1 


For 5 


4.6 ± 0.2 


8.5 


3.2 ±0.2 


-7.5 


4.1 ±0.3 



uncertainties due to the subtraction of the field star component 
from the spectral scans are very minor, even for Fornax 4. 

Binning of the CMDs. As can be seen in Fig. [4] our adopted 
binning is sufficient to well sample all regions of the CMDs. We 
verified that changes in the binning did not affect our results to 
any significant extent: increasing the number of cmd-bins fur- 
ther only changed the [Fe/H] values by about 0.01 dex. In fact, a 
somewhat smaller number of bins would probably be sufficient; 
iMc William & Bernsteinl (120081) used only 27 bins in their mod- 
elling of the integrated spectrum of the globular cluster 47 Tuc. 

Summary. From the above discussion, we conclude that the 
error sources discussed here all contribute by less than ~ 0.1 
dex to the uncertainties on our [Fe/H] measurements, and typi- 
cally only by a few times 0.01 dex or less to abundance ratios. 
Nevertheless, the true uncertainties may well be dominated by 
effects that we cannot easily quantify here, such as non-LTE line 
formati on, inclusion of sc attering in the synthetic spectral calcu- 
lations dCavrel et al.l2004l) . accuracy of damping parameters and 
the atomic oscillator strengths and the quality and completeness 
of the line lists in general, as well as the simplifications involved 
in the use of classical model atmospheres. Those limitations are 
not unique to our analysis, however, but affect most chemical 
abundance studies. 

5.4. Radial velocities 

The radial velocities of the clusters were determined with high 
accuracy as a byproduct of our analysis. The raw shifts required 
to match the data with the synthetic spectra were -72.4 km s _1 , 
-59.4 km s _1 and -72.7 km s _1 for Fornax 3, 4 and 5, respec- 
tively, with estimated uncertainties of ~0.2 km s _1 for Fornax 3 
and 5 and < 0.1 km s for Fornax 4. After applying heliocentric 
velocity corrections, this corresponds to radial vel ocities of 60.4 
km s" 1 , 47.2 km s _1 and 60.6 km s _1 . Previously, Dubatf Tet alJ 
(1992) have reported radial velocities of 58.5 km s -1 , 46.2 km 
s~ 2 and 62.1 km s _1 with errors of 1-2 km s , consistent with 
our measurements. 

The average radial velocities for Fornax 1, 2 and 3 measured 
bv lLetarte et~al] d2006l) are 59.3, 63.8 and 62.7 km s" 1 . However, 
these are based on measurements of three individual stars in each 
cluster and given the velocity dispersions of several km s _1 , these 
mean values must thus be considered uncertain by a few km s _1 . 

For comparison, the radial velocity of the Fornax dSph itself 
has been determined t o be 53.3 km s" 1 an d the velocity disper- 
sion o- = 11.1 km s _1 dMSsEiliDllQOl)- 

5.5. Velocity dispersions and dynamical masses 

As another byproduct of the analysis, we briefly discuss the dy- 
namics of the clusters. The cr smoot h//i ratios correspond to veloc- 
ity broadenings of ca smoo ±/A = 8.1 km s _1 , 5.5 km s _I and 6.3 
km s _1 . This may be compared with the (central) velocity dis- 



persions for Fornax 3, 4 and 5, which have been measured to 
be 8.8 km s _1 , 5. 1 km s~ l and 7.0 km s _1 , with uncertainties of 
about 1 km s 1 (iDubath etal.lll992l) . However, the instrumen- 
tal broadening needs to be taken into account. From a Gaussian 
fit to the sky emission lines at 5577 A and 5588/90 A we find 
a FWHM of 0.13 A, in excellent agreement with the resolving 
power /?=40000 expected for a 1" slit according to the UVES 
manual. This corresponds to crj nstr ~ 0.06 A or ~ 3 km s _1 . After 
subtracting this in quadrature from our measurements, we get 
cr v( j = 7.5 km s _1 , 4.6 km s _1 and 5.5 km s _1 for Fornax 3, 4 and 
5, respectively. Within the uncertainties, this ag rees fairly well 
with the measurements of Dubathetal. (1992). Our estimates 
tend to be systematically somewhat lower than theirs, but the er- 
rors on the Dubath measurements of the three clusters are prob- 
ably not entirely independent as they are based on the widths 
of cross-correlation functions from the CORAVEL instrument, 
corrected for instrumental resolution with a set of reference stars 
common to all three clus ters. By in tegrating the projected veloc- 
ity dispersion profiles of Kingj d 19661) models with concentration 
parameters in the range Wo = 5 — 7, we find that the central 
velocity dispersions are 5-10% higher than the values measured 
within the half-light radius, an effect that works in the right sense 
but is not sufficiently strong to fully explain the difference be- 
tween our measurements and those of Dubath et al. In any case, 
it is clear that the broadening of the spectral features in the GC 
spectra is dominated by the internal velocity dispersions for all 
three clusters. 

From the velocity dispersions and the half-light radii, we can 
derive dynamical masses, My u . We assume that 



where r;, is the three-dimensional half-mass radius, related to 
the projected two-dimensional half-light radius r e g as r/, a |r e ff 
(Spitzer 1987). This assumes that mass follows light, which may 
not necessarily be the case if the clusters are mass segregated. 
Further, cr^ denotes the global velocity dispersions. These are 
lower than our measured velocity dispersions <x vd , which only 
sample stars out to the half-light radii. We converted our velocity 
dispersion measurements within the rectangular regions covered 
by the UVE S observations to glo bal values by integ rating over 
iKingl d 19661) model profiles, as in lLarsen et alJ (120021) . 

In Table [8] we list the resulting global velocity dispersions, 
half-mass radii and virial masses as well as the V-band mass- 
to-light ratios (M/Ly). The quoted errors only take into account 
the uncertainties on the velocity dispersions, determined from 
the scatter around the red lines in F i g. El T he My magnitudes are 
based on the catalogue of Webbink ( 1985). Within t he errors, our 
masse s and M/Ly ratios agree well with those of Dubath et al. 
(Il992l) . In particular, we confirm that Fornax 4 has a signifi- 
cantly lower M/Ly ratio than the other, more metal -poor clus- 
ters. While this trend runs contrary to predictions by standard 
SSP models, it is similar to that seen for a much larger sample 
of GCs in M31 dStraderet all 120091 l201lh . Indeed, the M/L v 
ratio of Fornax 4 is very similar to those of old M3 1 globular 
clusters of similar metallicity. The M/Ly ratios of Fornax 3 and 
5 are very high, even by comparison with the mo st metal-poor 
cluste rs in M31 of which most have M/Ly < 3 dStrader et alJ 
1201 ll) . Dynamical effects will, in general, tend to lower the 
observed mass-to-light ratios relative to canonical SSP model 
predictions, both via the preferential loss of low-mass stars 
dKruijssen & Lamersll2008l) . and by leading to an underestimate 
of the true masses by dynamical measurements because the mea- 
sured velocity dispersions and structural parameters will be most 
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Table 9. Iron abundances, [Ca/Fe] ratios and radial velocities for 
the five GCs in Fornax from high-dispersion spectroscopy. 





[Fe/H] 


[Ca/Fe] 


v r (km s ') 


Source 


Fornax 1 


-2.5 ±0.1 


+0.15 ±0.04 


59 ± 1 


L06 


Fornax 2 


-2.1 ±0.1 


+0.20 ± 0.03 


64 ± 1 


L06 


Fornax 3 


-2.3 ±0.1 


+0.25 ± 0.08 


60.4 ± 0.2 


This work 


Fornax 4 


-1.4 ± 0.1 


+0.13 ±0.07 


47.2 ±0.1 


This work 


Fornax 5 


-2.1 ± 0.1 


+0.27 ± 0.09 


60.6 ± 0.2 


This work 



Notes. For [Fe/H] the overall errors are set to ±0.1 dex, including an 
estimate of systematic errors. For the remaining entries the errors are 
standard errors on the mean, based on the r.m.s. deviation of individual 
measurements. 



sensitive to the more massiv e stars, located preferentially near 
the centre (iFleck et alj 120061: lMiocchill2006l) . More exotic ex- 
planations include non-standard IMFs or even dark matter. Both 
top-heavy and bottom-heavy IMFs may in fact contribute to 
higher M/Ly ratios - the latter due to the increased number of 
faint low-mass stars , and the former due to a n increased number 
of dark remnants (iDabringhausen e t al. 2009). Binary stars could 
also lead to inflated velocity dispersions in integrated-light mea- 
surements. Common to all these potential explanations is that it 
is unclear why they should apply exclusively to the metal-poor 
Fornax GCs. 



6. Discussion 

6.1. Overall metallicities 

Combining the data from L06 and the present work, all five GCs 
in the Fornax dSph now have chemical abundance measurements 
from high-dispersion spectroscopy. In Table|9]we summarize the 
[Fe/H] and [Ca/Fe] values, where the latter may be taken as a 
proxy for the [a/Fe] ratio. For reference, we also include the 
radial velocity measurements, using our own measurements for 
Fornax 3, 4 and 5 and those of L06 for Fornax 1 and 2. From 
this table it is now clear that the metallicity of Fornax 4 is indeed 
much higher (by almost a full dex) than the average of the other 
clusters. Because our measurements do not rely on any interme- 
diate calibration steps of the metallicity scales, we believe this 
is a solid result. Moreover, any differences in the CMDs (such 
as the redder horizontal branch morphology of Fornax 4) are 
explicitly taken into account in our analysis. We note that the 
CMP of Fornax 4 is, in fact, fairly well fit by a iDotter et"ai] 
(120071) isochrone with [Fe/H] = -1.5 (Fig. g). This is partly 
he lped by the fact that we adopt a smaller E(V—I) value than that 
of Buonanno et al. ( 1999), making the RGB intrinsically redder 
and thus consistent with a higher metallicity. 

It is interesting to compare the metallicity distribution of 
the field stars in Fornax with the GC data. Field st ar metallic- 
ities have been measured by iBattaglia et alJ d2006l) from Ca n 
IR triplet spectroscopy. Comparing their metallicity distribution 
for the field stars (their Figure 19), we see that Fornax 1, 2, 3 
and 5 all belong at the extreme metal-poor tail of this distribu- 
tion. A more detailed analysis shows that a very large fraction, 
about 1/5-1/4, of the most metal-poor stars in Fornax are found 
within the four most metal-poor GCs dLarsen et al .1 120121) . The 
differences in chemistry are correlated with the spatial locations 
of the clusters, with Fornax 4 being located close to the centre 
of the Fornax dSph while th e other, more metal-p oor clusters are 
found further out, and in fact Strad er et al.l J2003) speculated that 
Fornax 4 may be the nucleus of the Fornax dSph. In this context 



it may be relevant that the radial velocity of Fornax 4 deviates by 
about 6 km s _I from the mean systemic velocity of the Fornax 
dSph. 

We may also compare the metallicities of the Forna x GCs 
with those in the Milky Way. The McMaster catalog dHarrisI 
1 19961) lists 1 1 clusters with [Fe/H] < -2, a fraction of only 7% 
of the Milky Way GC system. In striking contrast, 4 out of 5 
GCs in Fornax have [Fe/H] < -2. In fact, the more "metal- 
rich" cluster, Fornax 4, has a metallicity close to the peak of the 
metallicity distributio n of M ilky Way halo clusters, which is at 
[Fe/H] » -1.5 dZinnlll985l) . It thus appears unlikely that a ma- 
jority of the Milky Way GC system could have been accreted 
from Fornax-like dwarf galaxies. 

6.2. a-elements 

The abundances of different a-elements (O, Mg, Si, Ca, Ti) are 
generally found to correlate tightly with each other in field stars 
both in the Galaxy and in dSphs, with a tendency for the [Mg/Fe] 
ratio to be slightly highe r than that of the other a-elements in the 
dSphs dVenn et al.ll2004t) . It is therefore interesting that we find 
a significantly lower [Mg/Fe] ratio in the Fornax GCs compared 
to [Ca/Fe] and [Ti/Fe]. In fact, our [Mg/Fe] ratios are consis- 
tent with being roughly Solar in all three GCs. A similarly low 
[Mg/Fe] ratio compared to other g- elements has been found for 
the int egrated light of GCs in M31 dColucci et alj|2009l) and the 
LMC dColucci et alJl2012l) . Colucci et al. found this to be ac- 
companied by an elevated Al abundance. It is thus likely that we 
are detecting the signatures of the abundance anomalies of light 
elements that are well- known in Galactic GCs, s pecifically the 
Mg-Al anticorrelation dGratton etalJl200ll 120121) . It is perhaps 
not highly surprising that these abundance anomalies are present 
in extragalactic GCs, given their ubiquity in Galactic GCs, but 
the demonstration that they are detectable in integrated light is a 
very important step towards establishing whether these anoma- 
lies are unique to ancient GCs. Low-mas s open clusters in th e 
Milky Way do not display these patterns (Panc ino et al.l l2010). 
However, it is unclear whether this implies fundamentally dif- 
ferent origins for open and globular clusters, or simply reflects 
the large differences in the typical masses of objects belonging 
to either category that have been subject to detailed study. Any 
low-mass analogs of the surviving GCs have disrupted long ago, 
and are thus inaccessible to direct observations. A large frac- 
tion of the Galactic halo ma y indeed have formed originally i n 
now-dispersed clusters (e.g. lKruiissen & Portegies Zw art 2009), 
so th e near-absence of the a bundance anomalies in halo field 
stars dMartell & Grebe]||2010h is suggestive that these disrupted 
low-mass clusters were chemically more homogeneous inter- 
nally than the surviving GCs. The remaining question, then, is 
whether young analogues of globular clusters (at least in terms of 
their masses) that are forming in the present-day Universe share 
the same anomalous abundance patterns as their older cousins. 
Although young clusters with masses > 1O 5 M are rare or ab- 
sent injhe_I^o£<d_Gjwp 2 Jh^ 

ies dLarsen & Richtleril2000l iLarsenl 120091: IWhitmordl2003l) . If 
abundance anomalies are the norm in massive star clusters, a 
detailed analysis of the chemical composition of nearby young, 
massive clusters might reveal their presence also in such objects. 

The remaining ff-elements that we could measure in the 
Fornax GCs (Ca, Ti) are clearly super-solar, but slightly less 
so than the typical enhancement of ~ +0.4 dex seen i n Milky 
Way globular clusters and halo stars dMcWilliamll 1 9971) . In this 
respect, the Fornax GCs follow a trend seen also in the dSph 
field stars and in the Large Magellanic Cloud, which also tend 
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to be less a-enhanced than Milky Way stars of similar metal- 
lic!^ (iTolstov et alJl2003LIVenn et al.|l20rj4tlJohnson et al.ll2006t 
Kirb v et al.ll201 U IColucci et al.ll2.012b . For Fornax 3 and 5, both 
Ca and Ti are enhanced at the level of 0.25-0.30 dex, and there 
is a hint that Fornax 4 is even less a-enhanced, with [Ca/Fe] and 
[Ti/Fe] ratios of ~ 0.12 dex. This is interesting given the sug- 
gestion that Fo rnax 4 may have forme d a few Gyr later than the 
other clusters (Bu onanno et al.lll999l) . These authors also sug- 
gested that a non-standard chemical composition might be re- 
sponsible for differences between the metallicities inferred from 
the RGB slope and other methods. Although we do find a dif- 
ference, it is unclear whether a relatively modest difference of 
~ 0.1 dex between the [ff/Fe] ratios of Fornax 4 and the other 
clusters would be sufficient to explain this. 

6.3. Fe-peak elements 

We measured several Fe-peak elements, namely: Sc, Cr, Mn, Fe 
and Cu. In spite of the common label assigned to these elements, 
it is not clear that they are all formed in a single nucleosynthetic 
channel and they do not all scale in a simple way with the Fe 
abundance. Indeed, we found both Mn and Cu to be significantly 

depleted with respect to Fe, with [Mn/Fe] 0.3 to -0.5 and 

[Cu/Fe] ~ -0.6 to < -1, although it should be noted that our 
measurement of Cu is based on only a single line (5106 A). 
We found Solar or slightly sub-solar [Cr/Fe] values, while the 
[Sc/Fe] values are all consistent with Solar within the uncertain- 
ties. 

It is well established that Mn is depleted in metal-poor 
Milky Way halo stars and globula r clusters, including oj Cen, 

at the level of fMn/Fel 0.5 (lGrattonlll989t INissen et alJ 

2000; ISobeck et alll2006t iRomano et al.ll201 ll) . while at metal- 
licities [Fe/H] > -1 the [Mn/Fe] ratio gradually increases to 
Solar. A similar Mn depletion has been found in metal-poor 
field stars in dSph ga l axies, including Fornax, and in the LMC 
(IShetrone et all 120031: IColucci et al.1 120TI iNorth et al.l l20Tl . 
According to our measurements, the Fornax GCs are simi- 
lar to metal-poor stars in the Fornax dSph and elsewhere in 
this respect. The usual interpretation of the low [Mn/Fe] ra- 
tios in metal-poor stars is that yields are strongly metallicity- 
depende nt, although the actual formation site of Mn remains un- 
certain dTimmesetalJll995l: IShetrone etal1l2003l) . The higher 
Mn abundances d erived for Fornax 3 (L06) and for a few LMC 
globular clusters (IColucci et al.l 120 l"2b are then quite puzzling, 
and it would be desirable with more measurements in order to 
verify whether there is an offset between the [Mn/Fe] ratios in 
field- and globular cluster stars in some galaxies. We note that 
the Mn features that are often measured by other authors (such as 
the Mn i triplet near 6020 A) are mostly too weak in the Fornax 
GCs and our Mn abundances mainly rely on the less commonly 
used Mn i lines near 4770 A. These do not have hyperfine split- 
ting included in the line list we are using. However, the Kurucz 
website now includes line lists with hyperfine splitting for many 
more transitions, including these lines, and we have verified that 
the Mn abundances change very little (decreasing by 0.01-0.05 
dex, with the largest change for Fornax 4) if we use the new line 
lists. 

The behaviour of Cu is fairly similar to that of Mn. At 
[Fe/H] ~ -2, [Cu/Fe] ratios are typically in the range -0.9 
to -0.5 for both field and GC halo star s, with an increase 
in the rCu/Fe] ratio at higher [Fe/H] (ISneden et al.l 1 199 it 
Simmerer et al.l [20031) . In this respect, our measurements of 
Fornax 3 and 5 are similar to those of Milky Way stars with 



similar metallicity. However, Fornax 4 with its higher overall 
metallicity falls significantly below the trend found for Milky 
Way stars. At [Fe/H] = -1.4, the typical Cu abundance in Milky 

Way stars is [Cu/Fe] 0.5, albeit with a 0.1-0.2 dex scatter. 

Our lower [Cu/Fe] ratio of -0.8 for Fornax 4 mimics the patterns 
observed in the LMC and dSph field stars, which also tend to fall 
below the Milky Way [Cu/Fel vs. |"Fe/H] re lation (IShetrone et al.l 
120031: IShetronell2004l: (Johnson et al.ll2006l) . 

The Sc abundance generally ten ds to follow Fe, both i n 
metal-poor s tars in the Milky W ay dGratton & Snedenlll99ll). 
in the LMC dJohnson et al.ll2006|) and in dSphs (Shetrone et al.l 
120031) . However. INissen et al.l (l2000l) find [Sc/Fe] decreasing 
from [Sc/Fe] ~ +0.2 at [Fe/H] < -1 to [Sc/Fe] ~ at solar 
metallicity, and suggest that Sc (intermediate in atomic number 
between Ca and Ti) behaves like an a-element. Our data suggest 
no significant enhancement of the [Sc/Fe] ratio in the Fornax 
GCs. 

Concerning Cr, iJohnson et al.l (120061) find solar to slightly 
sub-solar [Cr/Fe] ratios for LMC GCs, and stars in the Milky 
Way appear to have [Cr/Fe] ratios very clos e to over a wide 
range of metallicity (Gratton & Sneden 1991). This is again sim- 
ilar to our results for the Fornax GCs. 



6.4. Neutron-capture elements 

We measured the abundances of the three neutron-capture ele- 
ments Ba, Y and Eu. Of these, Eu is considered a pure r-process 
element, while Ba and Y can also be produced via the s-process. 
As such, the abundance of Eu relative to the other elements is an 
indicator of the extent to which the nucleosynthesis was domi- 
nated by the r-process prior to the formation of the GCs. Since 
the site of the r-process is usually assumed to be massive stars, a 
strong r-process contribution is then taken as an indication that 
the stars in question formed on a short time scale. 

In the Milky Way halo, the neutron capture elements are 
characterized by a very large abundance spread relative to Fe, 
usually taken as evidence of inhomogeneous chemical evolu- 
tion. Notably, this spread is also seen internally in some GCs 
(Sect. 15.11 ). At metallicities [Fe/H] < -2, the Ba abundance 
spans a wide ran ge of -1 < [Ba/Fe ] < 1, with a few eve n more 
extreme outliers lMcWilliam|[l997l: iFrancois et al.ll2007l) . Such 
a spread is also seen in dSphs dVenn et al.ll2004l) . 

We found a rather large spread in the [Ba/Fe] ratios of the 
Fornax GCs, similar to that seen in the Galactic halo. Eu is also 
displaying a substantial spread, but our measurements of this el- 
ement are quite uncertain. Nevertheless, the Eu abundances in 
all three clusters appear strongly enhanced relative to both Fe 
and Ba with an average [Ba/Eu] 0.7. This [Ba/E u] ratio ap- 
proac hes that of the pure r-process ([Ba/Eu] = -0.8: iMc William! 
1997), indicating a dominant contribution to neutron capture 
synthesis by the r-process. This is again similar to the patterns 
observed in metal-poor stars in the Milky Way, a s well as in the 
LMC and Local Group d Sphs dVenn et al.ll2004t IJohnson et all 
120061: IColucci et ai1l2012l) . The [Ba/Eu] ratio of Fornax 4 is not 
significantly different from that of the other clusters - in fact, 
Fornax 4 formally has the lowest [Ba/Eu] ratio, and thus the 
strongest r-process signature, of the three. This may be some- 
what contrary to what we might expect if Fornax 4 were several 
Gyr younger than the other clusters, but the large uncertainties 
on the Eu abundances should be kept in mind here. They are 
based on measurements of two lines at 4205 A and 4436 A, both 
of which are rather complex blends. 
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The [Ba/Y] ratio is interesting as it measures the ratio of 
heavy vs. light neutron capture elements, and is less uncertain 
than the ratios involving Eu. We found a mean [Ba/Y] — 1-0.66, 
very simil ar to the mean [Ba/YJ ratio for field s tars in d Sphs 
and L MC (IShetrone et alj|2003l I Venn et alj|2004t IColucci et al.1 
1201 2|) . but significantly more po sitive than for metal-poor Milky 
Way stars dFrancois et all2007l) . 

7. Conclusions 

We have carried out a detailed analysis of metallicities and abun- 
dance patterns for the three globular clusters Fornax 3, 4 and 5 
in the Fornax dwarf spheroidal galaxy. Our analysis is based on 
the integrated light of the clusters, using drift-scan observations 
with the UVES spectrograph on the ESO Very Large Telescope, 
and a combination of Hubble Space Telescope observations of 
the colour-magnitude diagrams and theoretical models to char- 
acterize the stellar contents of the clusters. We have introduced 
and described the method we have developed for the analysis of 
such data, and carried out several tests to assess various system- 
atic and random errors. Our main results are as follows: 

- Combining our d ata with measurem ents of individual stars in 
Fornax 1 and 2 (iLetarte et al.l2006h . it is clear that the Fornax 
GCs fall into two distinct groups: Fornax 1, 2, 3 and 5 have 
very low metallicities in the range -2.5 < [Fe/H] < -2, well 
below the average metallicity of Milky Way halo GCs, while 
Fornax 4 has [Fe/H] 1.4. Taking into account both ran- 
dom and systematic errors, we estimate the uncertainties on 
the iron abundances to be ~ 0.1 dex. It is thus firmly estab- 
lished that Fornax 4 is significantly more metal-rich than the 
other four clusters. 

- Our integrated-light measurements of a wide range of ele- 
mental abundances in the Fornax GCs are mostly consistent 
with observations of individual metal-poor stars in dSphs, 
with the notable exception of [Mg/Fe]. Specifically, we find 
a more moderate a-enhancement (based on [Ca/Fe] and 
[Ti/Fe]) than typical of metal-poor stars in the Milky Way, 
especially for Fornax 4. We also find a strong r-process sig- 
nature, even in Fornax 4 despite the younger age suggested 
by some studies and the low [ayFe] ratio. 

- In all three clusters, the [Mg/Fe] ratio is significantly lower 
than [Ca/Fe] and [Ti/Fe]. This is likely a signature of the 
Mg-Al abundance anti-correlation that is well-known from 
some Milky Way globular clusters. While it is very encour- 
aging that such abundance anomalies may be detectable in 
integrated light, this also implies that [Mg/Fe] may be a poor 
proxy for the overall [a/Fe] ratio, at least for globular clus- 
ters. 

- We confirm that Fornax 4 has a significantly lower mass-to- 
light ratio than the more metal-poor clusters Fornax 3 and 
Fornax 5. This trend, while opposite to predictions by SSP 
models, is similar to that found for globular clusters in M3 1 
(IStrader et alJl2009ll20Tlh . 

Having the ability to carry out detailed abundance analysis 
from integrated light means that this type of analysis can now be 
extended to galaxies well beyond the Local Group. We believe 
that, for the time being, it is still risky to rely on purely theoreti- 
cal models for the underlying H-R diagrams of the clusters, par- 
ticularly at young ages where such models are still quite uncer- 
tain. However, colour-magnitude diagrams are available for sev- 
eral young ma ssive star clusters in galaxies within distances of 
several Mpc dLarsen et al.l 1201 ll) . allowing detailed integrated- 
light abundance analysis for such clusters to be put on a solid 



footing. In the slightly more distant future, the combination of 
high-resolution imaging and high-dispersion spectroscopy on 
30-40 m class telescopes will be a very powerful tool for de- 
tailed analysis of the chemical composition of star clusters in a 
significant volume of the Local Universe. 
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1.178 


1 


E 


0.207 


RGB/MS 


-1.764 


1.182 


1 


E 


0.213 


RGB/MS 


-1.719 


1.172 


1 


E 


0.218 


RGB/MS 


-1.692 


1.173 


1 


E 


0.223 


RGB/MS 


-1.639 


1.137 


1 


E 


0.227 


RGB/MS 


-1.637 


1.138 


1 


E 


0.232 


RGB/MS 


-1.537 


1.130 


1 


E 


0.236 


RGB/MS 


-1.403 


1.115 


1 


E 


0.240 


RGB/MS 


-1.366 


1.088 


1 


E 


0.243 


RGB/MS 


-1.317 


1.109 


1 


E 


0.247 


RGB/MS 


-1.304 


1.070 


1 


E 


0.250 


RGB/MS 


-1.284 


1.088 


5 


E 


0.267 


RGB/MS 


-1.096 


1.063 


5 


E 


0.281 


RGB/MS 


-1.000 


1.056 


5 


E 


0.294 


RGB/MS 


-0.945 


1.063 


5 


E 


0.306 


RGB/MS 


-0.780 


1.032 


5 


E 


0.317 


RGB/MS 


-0.639 


1.016 


5 


E 


0.326 


RGB/MS 


-0.444 


0.971 


5 


E 


0.334 


RGB/MS 


-0.304 


0.977 


5 


E 


0.341 


RGB/MS 


-0.227 


0.963 


5 


E 


0.347 


RGB/MS 


-0.179 


0.962 


5 


E 


0.353 


RGB/MS 


-0.138 


0.963 


5 


E 


0.359 


RGB/MS 


-0.097 


0.954 


5 


E 


0.365 


RGB/MS 


-0.077 


0.946 


5 


E 


0.370 


RGB/MS 


-0.021 


0.947 


5 


E 


0.376 


RGB/MS 


0.016 


0.923 


5 


E 


0.381 


RGB/MS 


0.118 


0.927 


5 


E 


0.385 


RGB/MS 


0.193 


0.916 


5 


E 


0.390 


RGB/MS 


0.252 


0.912 


5 


E 


0.394 


RGB/MS 


0.307 


0.923 


5 


E 


0.398 


RGB/MS 


0.327 


0.919 


5 


E 


0.401 


RGB/MS 


0.557 


0.900 


3.55e+01 


S 


0.423 


RGB/MS 


0.726 


0.886 


2.63e+01 


S 


0.437 


RGB/MS 


0.984 


0.879 


3.03e+01 


S 


0.450 


RGB/MS 


1.179 


0.848 


3.69e+01 


S 


0.463 


RGB/MS 


1.330 


0.834 


3.46e+01 


S 


0.473 


RGB/MS 


1.644 


0.822 


1.00e+02 


S 


0.496 


RGB/MS 



S. S. Larsen et al.: High-dispersion spectroscopy of the Fornax dSph GCs, Online Material p 2 
Table 1. continued. 



My 


V-I 


Weight 




f(L V ) 




1.963 


0.804 


1.17e+02 


S 


0.515 


RGB/MS 


2.195 


0.792 


8.26e+01 


S 


0.527 


RGB/MS 


2.327 


0.786 


8.17e+01 


S 


0.537 


RGB/MS 


2.486 


0.773 


8.97e+01 


S 


0.546 


RGB/MS 


2.632 


0.764 


9.89e+01 


S 


0.555 


RGB/MS 


2.765 


0.760 


7.52e+01 


S 


0.561 


RGB/MS 


2.975 


0.713 


2.62e+02 


S 


0.578 


RGB/MS 


3.213 


0.616 


4.22e+02 


S 


0.601 


RGB/MS 


3.420 


0.551 


4.90e+02 


S 


0.623 


RGB/MS 


3.550 


0.527 


4.07e+02 


S 


0.638 


RGB/MS 


3.677 


0.513 


4.81e+02 


S 


0.655 


RGB/MS 


3.779 


0.503 


4.65e+02 


S 


0.670 


RGB/MS 


3.883 


0.522 


5.92e+02 


s 


0.687 


RGB/MS 


3.982 


0.507 


5.86e+02 


s 


0.702 


RGB/MS 


4.069 


0.502 


4.68e+02 


s 


0.714 


RGB/MS 


4.142 


0.534 


4.82e+02 


s 


0.725 


RGB/MS 


4.220 


0.521 


5.32e+02 


s 


0.736 


RGB/MS 


4.294 


0.528 


5.63e+02 


s 


0.747 


RGB/MS 


4.379 


0.527 


5.65e+02 


s 


0.757 


RGB/MS 


4.456 


0.552 


5.14e+02 


s 


0.766 


RGB/MS 


4.665 


0.569 


1.57e+03 


s 


0.788 


ISO 


4.876 


0.592 


1.82e+03 


s 


0.809 


ISO 


5.090 


0.620 


2.24e+03 


s 


0.830 


ISO 


5.319 


0.651 


2.98e+03 


s 


0.853 


ISO 


5.572 


0.687 


3.95e+03 


s 


0.877 


ISO 


5.830 


0.727 


4.44e+03 


s 


0.899 


ISO 


6.049 


0.766 


3.87e+03 


s 


0.914 


ISO 


6.255 


0.802 


3.84e+03 


s 


0.926 


ISO 


6.486 


0.844 


4.73e+03 


s 


0.939 


ISO 


6.704 


0.886 


4.80e+03 


s 


0.949 


ISO 


6.944 


0.934 


5.65e+03 


s 


0.959 


ISO 


7.191 


0.984 


6.12e+03 


s 


0.967 


ISO 


7.448 


1.040 


6.69e+03 


s 


0.974 


ISO 


7.649 


1.085 


5.54e+03 


s 


0.979 


ISO 


7.851 


1.130 


6.13e+03 


s 


0.984 


ISO 


8.052 


1.173 


6.59e+03 


s 


0.988 


ISO 


8.264 


1.217 


7.75e+03 


s 


0.992 


ISO 


8.472 


1.257 


8.87e+03 


s 


0.995 


ISO 


8.735 


1.302 


1.39e+04 


s 


1.000 


ISO 



Table 2. Colour-magnitude diagram for Fornax 4 



M v 


V-I 


Weight 




f(Lv) 




0.583 


0.604 


14 


E 


0.009 


HB 


0.597 


0.668 


14 


E 


0.019 


HB 


0.572 


0.718 


14 


E 


0.028 


HB 


0.597 


0.749 


14 


E 


0.038 


HB 


0.505 


0.793 


14 


E 


0.048 


HB 


0.298 


0.859 




E 


0.051 


HB 


-2.617 


1.582 




E 


0.064 


RGB/MS 


-2.426 


1.450 




E 


0.075 


RGB/MS 


-2.414 


1.484 




E 


0.086 


RGB/MS 


-2.283 


1.408 




E 


0.095 


RGB/MS 


-2.269 


1.396 




E 


0.104 


RGB/MS 


-2.147 


1.351 




E 


0.113 


RGB/MS 


-2.129 


1.345 




E 


0.121 


RGB/MS 


-2.073 


1.288 




E 


0.129 


RGB/MS 


-2.064 


1.291 




E 


0.136 


RGB/MS 


-2.021 


1.310 




E 


0.144 


RGB/MS 


-1.937 


1.276 




E 


0.151 


RGB/MS 



S. S. Larsen et al.: High-dispersion spectroscopy of the Fornax dSph GCs, Online Material p 3 
Table 2. continued. 



M v 


V-I 


Weight 




f(L V ) 




-1.934 


1.246 


1 


E 


0.158 


RGB/MS 


-1.816 


1.265 




E 


0.164 


RGB/MS 


-1.817 


1.248 


1 


E 


0.170 


RGB/MS 


-1.723 


1.184 


1 


E 


0.175 


RGB/MS 


-1.693 


1.257 


1 


E 


0.181 


RGB/MS 


-1.657 


1.173 


1 


E 


0.186 


RGB/MS 


-1.627 


1.206 


1 


E 


0.191 


RGB/MS 


-1.584 


1.216 


1 


E 


0.196 


RGB/MS 


-1.568 


1.175 


1 


E 


0.201 


RGB/MS 


-1.461 


1.176 


5 


E 


0.223 


RGB/MS 


-1.210 


1.109 


5 


E 


0.241 


RGB/MS 


-1.117 


1.151 


5 


E 


0.257 


RGB/MS 


-1.045 


1.031 


5 


E 


0.272 


RGB/MS 


-0.862 


1.067 


5 


E 


0.285 


RGB/MS 


-0.758 


1.041 


5 


E 


0.297 


RGB/MS 


-0.584 


1.029 


5 


E 


0.307 


RGB/MS 


-0.450 


1.061 


5 


E 


0.315 


RGB/MS 


-0.364 


0.980 


5 


E 


0.323 


RGB/MS 


-0.184 


1.047 


5 


E 


0.330 


RGB/MS 


-0.086 


1.037 


5 


E 


0.336 


RGB/MS 


-0.015 


0.989 


5 


E 


0.342 


RGB/MS 


0.143 


1.028 


5 


E 


0.347 


RGB/MS 


0.166 


0.975 


5 


E 


0.352 


RGB/MS 


0.266 


0.993 


5 


E 


0.357 


RGB/MS 


0.298 


1.003 


5 


E 


0.361 


RGB/MS 


0.413 


0.993 


2.15e+01 


S 


0.378 


RGB/MS 


0.698 


0.944 


4.67e+01 


S 


0.406 


RGB/MS 


1.103 


0.929 


6.52e+01 


S 


0.434 


RGB/MS 


1.409 


0.883 


3.79e+01 


S 


0.446 


RGB/MS 


1.584 


0.875 


4.49e+01 


S 


0.458 


RGB/MS 


1.740 


0.875 


4.65e+01 


S 


0.469 


RGB/MS 


1.898 


0.866 


5.00e+01 


S 


0.479 


RGB/MS 


2.066 


0.873 


7.30e+01 


S 


0.491 


RIDGE 


2.263 


0.861 


8.87e+01 


S 


0.504 


RIDGE 


2.460 


0.848 


1.13e+02 


S 


0.517 


RIDGE 


2.677 


0.824 


1.30e+02 


S 


0.530 


RIDGE 


2.874 


0.780 


1.41e+02 


S 


0.542 


RIDGE 


3.070 


0.730 


1.64e+02 


S 


0.553 


RIDGE 


3.268 


0.662 


2.68e+02 


S 


0.568 


RIDGE 


3.465 


0.601 


5.59e+02 


S 


0.595 


RIDGE 


3.662 


0.557 


6.84e+02 


S 


0.622 


RIDGE 


3.879 


0.549 


8.89e+02 


s 


0.651 


RIDGE 


4.076 


0.553 


1.09e+03 


s 


0.680 


RIDGE 


4.277 


0.556 


1.33e+03 


s 


0.710 


ISO 


4.482 


0.571 


1.62e+03 


s 


0.740 


ISO 


4.728 


0.595 


2.26e+03 


s 


0.774 


ISO 


4.933 


0.617 


2.14e+03 


s 


0.800 


ISO 


5.147 


0.644 


2.54e+03 


s 


0.825 


ISO 


5.369 


0.674 


2.98e+03 


s 


0.850 


ISO 


5.593 


0.705 


3.36e+03 


s 


0.872 


ISO 


5.819 


0.737 


3.72e+03 


s 


0.893 


ISO 


6.040 


0.772 


3.93e+03 


s 


0.910 


ISO 


6.264 


0.811 


4.22e+03 


s 


0.925 


ISO 


6.511 


0.856 


4.94e+03 


s 


0.939 


ISO 


6. /30 


a oao 


A C A ~ i 

4.54e+03 


c 


A ACA 

0.950 


TCA 

1SU 


6.973 


0.949 


5.17e+03 


s 


0.960 


ISO 


7.238 


1.011 


5.82e+03 


s 


0.968 


ISO 


7.450 


1.063 


4.79e+03 


s 


0.974 


ISO 


7.650 


1.116 


4.66e+03 


s 


0.979 


ISO 


7.891 


1.180 


6.02e+03 


s 


0.983 


ISO 


8.138 


1.245 


6.89e+03 


s 


0.988 


ISO 
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Table 2. continued. 



M v V-I Weight f(L v ) 

8.418 1.314 8.87e+03 S 0.992 ISO 

8.708 1.378 1.12e+04 S 0.997 ISO 

8.954 1.426 1.15e+04 S 1.000 ISO 



Table 3. Colour-magnitude diagram for Fornax 5 



M v 


V-I 


Weight 




f(Lv) 




1.064 


-0.113 


14 


E 


0.011 


HB 


0.813 


-0.054 


14 


E 


0.025 


HB 


0.657 


0.017 


14 


E 


0.041 


HB 


0.600 


0.072 


14 


E 


0.058 


HB 


0.506 


0.145 


14 


E 


0.077 


HB 


0.478 


0.286 


14 


E 


0.096 


HB 


0.494 


0.391 


14 


E 


0.114 


HB 


0.442 


0.545 


14 


E 


0.134 


HB 


0.063 


0.809 


14 


E 


0.162 


HB 


-0.620 


0.944 


14 


E 


0.214 


HB 


-2.554 


1.322 


1 


E 


0.236 


RGB/MS 


-2.509 


1.291 


1 


E 


0.258 


RGB/MS 


-2.233 


1.250 


1 


E 


0.274 


RGB/MS 


-2.198 


1.225 


1 


E 


0.290 


RGB/MS 


-2.138 


1.215 


1 


E 


0.305 


RGB/MS 


-2.023 


1.185 


1 


E 


0.319 


RGB/MS 


-1.821 


1.163 


1 


E 


0.330 


RGB/MS 


-1.480 


1.117 


1 


E 


0.338 


RGB/MS 


-1.393 


1.113 


1 


E 


0.346 


RGB/MS 


-1.384 


1.088 


1 


E 


0.353 


RGB/MS 


-1.336 


1.075 


1 


E 


0.361 


RGB/MS 


-1.323 


1.061 


1 


E 


0.368 


RGB/MS 


-1.283 


1.060 


1 


E 


0.375 


RGB/MS 


-1.144 


1.059 


1 


E 


0.381 


RGB/MS 


-1.108 


1.042 


1 


E 


0.387 


RGB/MS 


-1.087 


1.044 


1 


E 


0.392 


RGB/MS 


-1.028 


1.030 


1 


E 


0.398 


RGB/MS 


-0.994 


1.052 


1 


E 


0.403 


RGB/MS 


-0.981 


1.044 


1 


E 


0.408 


RGB/MS 


-0.910 


1.026 


1 


E 


0.413 


RGB/MS 


-0.736 


1.011 


5 


E 


0.434 


RGB/MS 


-0.568 


0.973 


5 


E 


0.452 


RGB/MS 


-0.010 


0.941 


5 


E 


0.462 


RGB/MS 


0.043 


0.928 


5 


E 


0.472 


RGB/MS 


0.071 


0.933 


5 


E 


0.482 


RGB/MS 


0.164 


0.920 


5 


E 


0.491 


RGB/MS 


0.234 


0.915 


5 


E 


0.500 


RGB/MS 


0.297 


0.920 


5 


E 


0.508 


RGB/MS 


0.381 


0.908 


5 


E 


0.515 


RGB/MS 


0.479 


0.897 


5 


E 


0.522 


RGB/MS 


0.580 


0.876 


5 


E 


0.528 


RGB/MS 


0.696 


0.897 


5 


E 


0.534 


RGB/MS 


0.776 


0.874 


5 


E 


0.539 


RGB/MS 


0.878 


0.875 


5 


E 


0.544 


RGB/MS 


0.937 


0.854 


5 


E 


0.548 


RGB/MS 


0.966 


0.879 


5 


E 


0.552 


RGB/MS 


1.014 


0.869 


5 


E 


0.557 


RGB/MS 


1.067 


0.837 


5 


E 


0.561 


RGB/MS 


1.144 


0.845 


5 


E 


0.564 


RGB/MS 


1.206 


0.844 


5 


E 


0.568 


RGB/MS 


1.347 


0.835 


2.16e+01 


S 


0.581 


RGB/MS 


1.567 


0.817 


3.26e+01 


S 


0.597 


RGB/MS 


2.165 


0.791 


9.80e+01 


S 


0.625 


RGB/MS 
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Table 3. continued. 



My 


V-I 


Weight 




f(L V ) 




2.515 


0.764 


7.54e+01 


S 


0.641 


RGB/MS 


2.791 


0.737 


8.34e+01 


S 


0.654 


RGB/MS 


2.980 


0.719 


5.37e+01 


S 


0.662 


RGB/MS 


3.150 


0.665 


8.47e+01 


S 


0.671 


RGB/MS 


3.264 


0.614 


6.74e+01 


S 


0.679 


RGB/MS 


3.336 


0.609 


5.82e+01 


S 


0.684 


RGB/MS 


3.397 


0.565 


6.49e+01 


S 


0.690 


RGB/MS 


3.463 


0.543 


7.51e+01 


S 


0.697 


RGB/MS 


3.529 


0.531 


7.46e+01 


S 


0.703 


RGB/MS 


3.658 


0.508 


2.99e+02 


S 


0.724 


RGB/MS 


3.839 


0.519 


3.78e+02 


S 


0.748 


RGB/MS 


4.014 


0.500 


3.96e+02 


S 


0.768 


RGB/MS 


4.167 


0.513 


3.62e+02 


S 


0.785 


RGB/MS 


4.291 


0.528 


3.52e+02 


S 


0.799 


RGB/MS 


4.422 


0.539 


3.13e+02 


S 


0.810 


RGB/MS 


4.528 


0.567 


2.82e+02 


S 


0.820 


RGB/MS 


4.631 


0.567 


2.78e+02 


S 


0.828 


RGB/MS 


4.723 


0.572 


2.75e+02 


s 


0.835 


RGB/MS 


4.813 


0.587 


2.64e+02 


s 


0.842 


RGB/MS 


4.904 


0.603 


3.40e+02 


s 


0.850 


RGB/MS 


4.984 


0.609 


3.24e+02 


s 


0.857 


RGB/MS 


5.202 


0.635 


9.85e+02 


s 


0.874 


ISO 


5.443 


0.668 


1.35e+03 


s 


0.893 


ISO 


5.703 


0.707 


1.67e+03 


s 


0.911 


ISO 


5.951 


0.748 


1.69e+03 


s 


0.926 


ISO 


6.170 


0.787 


1.55e+03 


s 


0.937 


ISO 


6.389 


0.826 


1.67e+03 


s 


0.947 


ISO 


6.591 


0.864 


1.67e+03 


s 


0.955 


ISO 


6.822 


0.909 


2.05e+03 


s 


0.963 


ISO 


7.067 


0.959 


2.30e+03 


s 


0.971 


ISO 




1 .viz 




c 


yj.y 1 1 


Ten 

lOVJ 


7.582 


1.070 


2.77e+03 


s 


0.982 


ISO 


7.851 


1.130 


3.17e+03 


s 


0.987 


ISO 


8.052 


1.173 


2.56e+03 


s 


0.990 


ISO 


8.264 


1.217 


3.02e+03 


s 


0.993 


ISO 


8.472 


1.257 


3.45e+03 


s 


0.996 


ISO 


8.735 


1.302 


5.39e+03 


s 


1.000 


ISO 



